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Abstract 

Scalar and tensor interactions were once competitors to the now well-established 
V — A structure of the Standard Model weak interactions. We revisit these inter- 
actions and survey constraints from low-energy probes (neutron, nuclear, and pion 
decays) as well as collider searches. Currently, the most stringent limit on scalar 
and tensor interactions arise from 0+ — > 0+ nuclear decays and the radiative pion 
decay tt — > eu^, respectively. For the future, we find that upcoming neutron beta 
decay and LHC measurements will compete in setting the most stringent bounds. 
For neutron beta decay, we demonstrate the importance of lattice computations of 
the neutron-to-proton matrix elements to setting limits on these interactions, and 
provide the first lattice estimate of the scalar charge and a new average of existing 
results for the tensor charge. Data taken at the LHC is currently probing these 
interactions at the 10~^ level (relative to the standard weak interactions), with 
the potential to reach the < 10^^ level. We show that, with some theoretical as- 
sumptions, the discovery of a charged spin-0 resonance decaying to an electron and 
missing energy implies a lower limit on the strength of scalar interactions probed 
at low energy. 



1 Introduction 



Nuclear and neutron beta decays have historically played a central role in determining the 
V—A structure of weak interactions and in shaping what we now call the Standard Model 
(SM) [HE]- Nowadays, precision measurements of low-energy processes such as neutron 
decay can be used to probe the existence of non-SM interactions, such as novel scalar 
and tensor structures. Considerable experimental efforts using both cold and ultracold 
neutrons are underway worldwide, with the aim to improve the precision of various 
neutron decay observables [31 E] : lifetime [5], El [71 El E] , beta asymmetry A [TOl [TH [121 [13] 
neutrino asymmetry B p3l [12], electron-neutrino correlation a [IHl [161 E]; Fierz 
interference term b [T5l|T8]. In some of the asymmetry measurements there are prospects 
to reach experimental sensitivities between 10^'^ and 10~^; this makes these observables 
very interesting probes of new physics effects originating at the TeV scale that have 
expected size (v/Absm)^, where v = {2\/2Gf)~^^'^ ~ 174 GeV and Absm denotes the 
mass scale where physics beyond the Standard Model (BSM) appears. 

The overall goal of this work is to assess the discovery potential and discriminating 
power of planned precision beta-decay measurements with cold and ultracold neutrons. 
In particular we wish to study the sensitivity of neutron decay to new physics in the 
context of and in competition with: (i) other low-energy precision measurements in 
nuclear beta decays and pion decays; and (ii) high-energy collider searches (Tevatron, 
LHC). In order to achieve our goal, we work within an effective field theory (EFT) setup, 
in which the dynamical effects of new heavy BSM degrees of freedom are parameterized 
by local operators of dimension higher than four built with SM fields. In the absence 
of a clear new-physics signal from collider searches, we find this way of proceeding the 
most attractive and general: all specific model analyses of beta decays (see Ref. [19] for a 
discussion within supersymmetry) can be cast in the EFT language and the constraints 
on effective operators that we will derive can be readily converted into constraints on 
the parameters of any SM extension. 

Among various BSM contributions we identify new scalar and tensor operators involv- 
ing left-handed neutrinos as the most promising to probe with neutron decay, because 
they interfere with the SM amphtude and thus contribute at linear order to decay pa- 
rameters. Motivated by this, in the unified EFT framework we present a comprehensive 
analysis of constraints on such scalar and tensor BSM interactions from a broad range 
of low-energy probes (neutron decay, nuclear decays, pion decays) as well as collider 
searches]^ To our knowledge such an analysis is missing in the literature, despite being 
essential to judging the relative merits of various low-energy experiments. 

Extracting bounds on short- distance scalar and tensor couplings from neutron and 
nuclear beta decays requires knowledge of the nucleon scalar and tensor form factors 
at zero momentum transfer, denoted here by gs^r- In previous beta-decay studies, gs 
and gT have been assumed to be 0(1) based on quark-model estimates (see, for example, 
Ref. 120]). The importance of the hadronic form factors can be appreciated by considering 
the extreme case in which gs,T <^ 1, which would dilute the sensitivity of beta decays to 
new physics. Concerning the hadronic form factors, the main results of this work are: 

-•^The EFT analysis of collider searches is valid as long as the particles that mediate the new interac- 
tions are above threshold for production at colliders. 
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• We provide the first lattice-QCD estimate of gs and a new average of existing qt 
results. Current lattice uncertainties are at the level of 50% for gs and 35% for qt- 
This already enables much improved phenomenology (see for example Fig. |9]) . 

• We show that a precision of 20% in gs will be needed to take full advantage of 10~^- 
level neutron-decay measurements. We identify and discuss the key systematic 
effects that need to be brought under control in order to achieve Sgs/gs ~ 20%. 

Besides the new estimates of gs and gT with lattice QCD (LQCD), the main new 
findings of our analysis can be summarized as follows: 

• Currently, the most stringent constraints on the scalar and tensor effective cou- 
plings (denoted by and er) arise from low-energy probes, is constrained by 
O"*" — )■ O"*" nuclear beta decays, while is constrained by the Dalitz-plot analysis 
of the radiative pion decay vr — )■ ew-f. There are also potentially very strong con- 
straints on es^T from the ratio of tt — )■ ez/ to vr — )■ /iz/ decay rates. This constraint 
arises from operator mixing: once a scalar or tensor interaction is generated by new 
physics, SM radiative corrections will generate an effective pseudoscalar operator 
that mediates the helicity-suppressed mode vr — )■ eu. If the flavor structure of the 
SM extension is known, this constraint could be the strongest. 

• Future neutron-decay measurements of the Fierz interference term b and the ana- 
logue term by in the neutrino asymmetry B can greatly improve existing constraints 
on tensor interactions: precision levels 6b, 56^ ~ 10^^ would provide a four-fold or 
higher improvement in the bound (depending on the sign of e^), as shown in Figs. |2] 
and [oj On the other hand, 6b, 6bi, ~ 10~^-level measurements would improve cur- 
rent bounds on er by one order of magnitude and current bounds on £5 by a factor 
of two (see Fig. Is] 



Current collider bounds from the LHC are not yet competitive with low-energy 



constraints (see Fig. 10). Folding in the current uncertainty on gs, the LHC bounds 
on es and er are weaker by a factor of about 4 and 3, respectively, than those 
obtained from nuclear decays and vr — t- ev^. 

Future LHC results, based on higher center-of-mass energy and higher integrated 
luminosity, would definitely improve on current low-energy bounds on es^r-, and 
would compete with improved low-energy constraints based on 6b, 6by ~ 10~^ in 



future neutron-decay measurements (see Fig. 11). 



Finally, we have explored the possibility that a mediator of new scalar interactions 
can be produced at the LHC. In this case, the EFT approach breaks down at 
collider energies and we have derived a general correlation between production 
cross-section for a scalar resonance at colliders and new-physics signal in neutron 
decay. This correlation links the discovery of a scalar resonance in — )■ ez/ + X 
at the LHC with a lower bound (i.e. guaranteed signal) on e^. This is illustrated 



in Figs. 12, 13, and 14 
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The paper is organized as follows. In Section[2]we present the effective theory descrip- 
tion of low-energy charged- current processes and briefly discuss how the coefficients may 
be constrained. In Section |3] we explain our notation for the matrix elements required to 
describe the neutron beta decay and discuss how this decay constrains the parameters 
in the effective field theory. In Section |4| we discuss the low-energy phenomenological 
constraints on chirality- violating scalar and tensor operators in the effective Lagrangian. 
Section [5] discusses current and planned lattice analyses for the matrix elements of the 
quark bilinear structures uTd between neutron and proton states, with special emphasis 
on the scalar and tensor structures. We provide the first estimate of gs from lattice QCD 
and a new average of existing calculations of qt- In Section [6] we summarize the impact 
of lattice estimates of gs,T on the phenomenology of scalar and tensor BSM interactions. 
In Section [7} we present the constraints on the short- distance couplings obtained from 
an analysis of high-energy scattering experiments and discuss the improvement expected 
in the next few years. We present our concluding remarks in Section [8j Two appendices 
provide details of the operators contributing to charged-current processes and of the 
neutron-decay differential decay distribution. 

2 Effective theory description of low-energy charged- 
current processes 

Following Ref . [21] , we describe new physics contributions to low-energy charged-current 
(CC) processes in a model-independent effective-theory setup, paying special attention 
to neutron-decay observables and their interplay with other low-energy and collider mea- 
surements. 

We parameterize the effect of new degrees of freedom and interactions beyond the 
SM via a series of higher-dimensional operators constructed with low-energy SM fields, 
assuming the existence of a mass gap between the SM and its ultraviolet completion. If 
the SM extension is weakly coupled, the resulting TeV-scale effective Lagrangian linearly 
reahzes the electroweak (EW) symmetry SU{2)l x U{1)y and contains a SM-like Higgs 
doublet [22]. We also assume that potential right-handed neutrino fields (sterile with 
respect to the SM gauge group) are heavy compared to the weak scale and therefore have 
been integrated out of the low-energy effective theory. This method is quite general and 
allows us to study the implications of precision measurements on a large class of models. 

In our analysis we truncate the expansion of the effective Lagrangian to the low- 
est non-trivial order, given by dimension-six operators. The contribution from the 
dimension-six operators to physical amplitudes involves terms proportional to f^/A|gjy[ 
and i?^/A|g(^, where v = {ip^) ~ 174 GeV is the vacuum expectation value (VEV) of the 
Higgs field and E is the characteristic energy scale of a given process. We will work to 
linear order in these ratios of scales. 

2.1 Effective Lagrangian 

In Ref. [21] a minimal basis of SU (2) x f/(l) invariant dimension-six operators contribut- 
ing to low-energy charged-current processes was identified (see Appendix |A] for details). 
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Denoting with Aj the effective dimensionful couphng associated with the operator Oj, 
we can write the effective Lagrangian as 



11 v"^ 
£(cff) ^ + ^ _ — y + _ ^ Oi , with = ^ , (1) 

i * i * 

where in the last step we have set the correct dimensions by the Higgs VEV v and defined 
the dimensionless new-physics couphngs dj, which are O(10~^) for Aj ~ TeV. 

In this framework one can derive the low-scale 0(1 GeV) effective Lagrangian for 
semi-leptonic transitions. It receives contributions from both VT-exchange diagrams 
(with modified VT-fermion couplings) and the four-fermion operators Oqde, Oiq, O* 
defined in Appendix |A| This matching procedure leads to [21] 



Iq 



1 + [vL]mj ) hlti^iL UiYd^L + [vnlmj h'Jt^i'eL UnYd' 



+ [sLluij ^R^iL U'j^d^ + [sR]uij ^rVrl u^d^R 



+ [tLlUij ^ROf.yl'iL URa'"'d{ 



+ h.c. . (2) 



where we use a^'^ = i['j^,'y'-']/2. The SM effective Lagrangian corresponds to vl = vr = 
sl = sr = ti = 0. The effective couplings vl, vr, sl, sr, ti ~ I are functions of the 
coupling OLi of SU{2) x U{1) invariant weak-scale operators. While their explicit expres- 
sions can be found in Appendix |X| here we simply point out two important features: 

• vl involves a linear combination of three weak-scale effective couplings: a quark- 
gauge boson vertex correction, a lepton-gauge boson vertex correction, and a four- 
fermion operator coupling left-handed quarks and leptons (same chirality structure 
as the SM). An important consequence is that by SU{2) x f/(l) gauge invariance, 
Vl is related to fermion-antifermion vertex corrections and neutral-current four- 
fermion vertices. 

• Vr and SL,SR,tL are in one-to-one correspondence with weak-scale effective cou- 
plings. Vr describes a right-handed charged- current quark coupling, while sl, sr, II 
correspond to scalar and tensor four-quark operators. Again, SU{2) gauge invari- 
ance implies that these couplings mediate not only charged-currrent processes but 
also processes such as ee uu, dd, with scalar or tensor Dirac structure. 

In what follows, we will work in the limit in which the effective non-standard couplings 
vl,r, sl,r, and ti are real and we will focus only on CP-even observables (for a discussion 
of CP-odd observables refer to Ref. [SO])- To simplify the notation, we will omit flavor 
indices, e.g. [vL]eeud vl- In addition, we will use the tree-level definition of the Fermi 
constant (7^/(8M^) = G'P /\/2. Working to linear order in the non-standard couplings, 
and focusing on the ij = ud component, the semi-leptonic effective Lagrangian can be 



4 



written in the following useful form: 



= -^J^ \ l + eL + en 



+ ^{l - 75)z/£ ■ u 



[I - 75)z/^ ■ u 



(1 - 2efi)7^75 



d 



es - ep75 



d + er iUf,u{l - ■ u(t'"'{1 - 75)^ 



+ I1.C., (3) 



where we have defined the effective scalar, pseudoscalar, and tensor couplings as follows: 



(^L,R = vl,r (^s = sl + sr ep = sl- sr ex = ti 



(4) 



While the physical amplitudes are renormalization scale and scheme independent, the 
individual effective couplings and hadronic matrix elements can display a strong scale 
dependence. Throughout the paper, we will quote estimates and bounds for the ej at 
the renormalization scale /x = 2 GeV in the MS scheme, unless otherwise specified. 

The Lagrangian ^ mediates all low-energy charged-current weak processes involving 
up and down quarks. For a recent analysis of fiavor-dependent constraints, see Ref. [23j . 
In some of the charged-current processes involving first-generation quarks the theoretical 
and experimental precision has reached or will reach in the near future a level that allows 
stringent bounds on the new-physics effective couplings. In this work we are interested 
in assessing the sensitivity of neutron decay to new physics in the context of (i) other 
low-energy constraints from nuclear beta decays and pion decays; and (ii) constraints 
from high-energy colliders (LEP, Tevatron, LHC). To set the stage for the discussion, we 
summarize the observables that give us access to the couplings appearing in Eq. ([s]) (we 
will come back in detail to these in following sections): 

• The combination (e/, + en) affects the overall normalization of the effective Fermi 
constant. This is phenomenologically accessible through quark-lepton universal- 
ity tests (precise determination of Vud from O''" — )■ 0+ nuclear decays under the 
assumption that Gp = G^, where is the Fermi constant extracted from muon 
decay). An extensive analysis of the constraints on {ep + ^R) from universality tests 
and precision electroweak observables from the Z-pole was performed in Ref. |21j . 
within BSM scenarios with minimal flavor violation. In this context it was shown 
that constraints from low-energy are at the same level or stronger (depending on 
the operator) than from Z-pole observables and e^e~ — )■ qq cross-section measure- 
ments at LEP. 

• The right-handed coupling ep affects the relative normalization of the axial and 
vector currents. In neutron decay ep can be reabsorbed in a redeflnition of the axial 
coupling and experiments are only sensitive to the combination (1 — 2eR)gA/gv 
{gv and gA are the vector and axial form factors at zero momentum transfer, to 
be precisely deflned below). Disentangling ep requires precision measurements of 
(1 — 2eR)gA/gv and precision calculations of gA/gv in LQCD. 

• The effective pseudoscalar combination ep = sp — sr contributes to leptonic 
decays of the pion. It is strongly constrained by the helicity-suppressed ratio 
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= r(7r ez/[7])/r(7r ^ /iz/[7]). Moreover, as discussed in Refs. pij 1251 126]. 
the low-energy coupling ep receives contributions proportional to es^r through elec- 
troweak radiative corrections. We will discuss the resulting constraints on es^p^r in 
Section STl 

• Both the scalar combination es = sl + sr and the tensor coupling ex = ti con- 
tribute at linear order to the Fierz interference terms in beta decays of neutrons 
and nuclei, and the neutrino-asymmetry correlation coefficient B in polarized neu- 
tron and nuclear decay (see Appendix [B] for notation). Because of the peculiar way 
in which the Fierz interference term appears in many asymmetry measurements, 
bounds on es and can also be obtained by observation of the beta-asymmetry 
correlation coefficient A, electron-neutrino correlation a, and positron polariza- 
tion measurements in various nuclear beta decays. Finally, the tensor coupling 
ex can also be constrained through Dalitz-plot studies of the radiative pion decay 
TT — )■ ez/7. 

• All of the above operators can provide signatures at colliders. Currently there are 
no competitive collider bounds on the chirality-flipping scalar and tensor couplings 
^s,p,T, because their interference with the SM amplitude carries factors of mf/Ef 
(where m/ is a light fermion mass, / G {e, u, d}), which at collider energies strongly 
suppresses the whole effect. So we immediately see that low-energy physics pro- 
vides a unique opportunity to probe these couplings, to which collider searches are 
sensitive only quadratically (i.e. via non-interference terms). We will derive in 
Section [7] the current bounds on es,T from searches at the LHC, and we will show 
that with higher center-of-mass energy and integrated luminosity they will become 
competitive with low-energy searches. 

Next, we review the analysis of neutron decay in the SM and beyond within the EFT 
framework described above. 

3 Neutron f3 decay 

The amplitude for neutron decay n{pn) — )• v{Vp)^ {'Pe)j^e{pu) mediated by the effective 
Lagrangian ([s]) involves in principle the matrix elements between the neutron and proton 
of all possible quark bilinears. These can be parameterized in terms of Lorentz-invariant 
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form factors as follows |27j : 



n{pn)) 



9T(v){q^) 

2Mm 



c'tiuq" + 



Uq') 



2Mn 



q^l UniPn) 



(5a) 



n{Pn)) 



Up{Pp) gAiq^hf,+ 



gT{A){q^) 
2Mn 



gp{q^) 



2Mn 



(5b) 



{p{pp)\ud 
{p{Pp)\ ui5d 
{p{Pp) \ ua^yd 



n{Pn)) 
n{Pn)) 

n{Pn)) 




+ gT\q^) (qt^Py - q^Pt^) + gPiq^) [i4iv - ivpy) Mn(pn) 



(5e) 



where the proton and neutron spinor amplitudes, P = pn + Pp, q = Pn — Pp is 

the momentum transfer, and M^v = M„ = Mp denotes a common nucleon mass g Note 
that all the above spinor contractions are 0(1), except for Up^y^Un which is 0{q/M]\i). 
Moreover, as discussed below, second-class current contributions gs and gT{A) affect the 
amplitude at levels below the expected experimental sensitivities. 

Our goal here is to identify TeV-induced new physics contaminations to the amplitude 
of typical size ep^s,T ~ {v/-^bsmY ~ 10~^. The effect we are after is of the same size 
as recoil corrections q/M^ ~ 10^^ as well as radiative corrections a/vr. So in our 
analysis we perform a simultaneous expansion in new physics contributions, recoil, and 
radiative corrections keeping terms up to first order and neglecting higher-order terms, 
as they are smaller than the current and planned experimental sensitivity. In light of this 
simultaneous expansion in ep^s,T, q/M^, and a/n, we now discuss contributions from all 
quark-bilinear operators: 

• Vector current: The form factor gv{0) contributes at 0(1) to the amplitude and 
^T(v)(0) contributes at first order in q/M^. Also, up to isospin-breaking corrections 
of order {Mn — Mp)/Mp^ ~ q/M]\f, the weak magnetism form factor ^T(y)(0) can be 
related to the difference of proton and neutron magnetic moments, that are well 
known. On the other hand, the induced-scalar form factor gs{q^) vanishes in the 
isospin limit [27], so it is of order (M„ — Mp)/MN ~ q/M^. Since it multiplies 
one power of q^/Mj^, its contribution to the amplitude is effectively second order 
in the recoil expansion, so we drop it. 

• Axial current: From the axial current only (7^(0) contributes up to first order. 
The induced-tensor form factor gT{A){q^) vanishes in the isospin limit [27j, and 
since it multiplies one power of q^/M^ its contribution to the amplitude is of 

^In the case of vector and axial bilinears, the induced tensor term proportional to (T^v<f can be 
traded for an independent "scalar" form factor proportional to P^. Here we choose to follow the 
parameterization of Ref. [27] . 



7 



second order in qjMj^^ so we drop it. Similarly, the contribution associated with 
the induced-pseudoscalar form factor is quadratic in our counting, because the 
pseudoscalar bilinear is itself of order g/M^r, and it comes with an explicit q/Mjq 
suppression, so we neglect itj^ 

Pseudoscalar bilinear: The pseudoscalar bilinear Up'y^Un is itself of order q/M^- 
Since it necessarily multiplies a new-physics effective coupling ep (there is no pseu- 
doscalar coupling in the SM), this term is also of second order in our expansion, 
and we drop it. 

Scalar and tensor bilinears: These bilinears enter into the analysis multiplied 
by new-physics effective couplings es,T- So we need the matrix elements to zeroth 
order in the recoil expansion, which leaves us with gs{0) and gri^)- gT''^''^\(f) are 
all multiplied by one power of q and gj.' vanishes in the isospin limit 



In summary, to the order we are working, the amplitudes depend only on gi = gi{Q) 
{i G {V, A, S,T}) and ^r(y)(0). Up to second-order corrections in isospin breaking, one 
has gv = ^ [291 [30] . For notational convenience, it is also useful to define the ratio 
of the axial to vector form factors a.s X = gA/gv- As noted earlier, in presence of 
non-standard right-handed interactions the axial form factor is always multiplied by the 
correction factor (1 — 2e/j), so that the neutron-decay amplitude is actually a function 
of A = A(l -2e/j). 

Finally, in order to make contact with the existing standard references on neutron 
and nuclear beta-decay phenomenology [311 1321 E], let us note here that Eq. (|5) can 
be viewed as the matching conditions from our quark-level effective theory Eq. (31) to 
a nucleon-level effective theory, such as the one originally written down by Lee and 
Yang [5T]. The Lee- Yang effective couplings Cj, C'^ {i G {V, A, S,T}) can be expressed 
in terms of our parameters as 

Q = ^VudCi (6a) 
Cv = gvH + tL + tR) (6b) 

Ca = -gA (1 + Cl - Cij) (6c) 

Cs = gs^s (6d) 
Ct = 4:gTeT , (6e) 

with Cj' = Ci, since we only have left-handed neutrinos in our low-energy effective theory. 
Operators involving right-handed neutrinos do not interfere with the SM amplitude 
and therefore contribute at second order to all observables. An analysis involving such 
operators will be presented elsewhere [33]. Finally, notice that Ref. [20] defines the 
couplings Ca, C(/5t with an overall minus sign compared to ours. 



■^This effect is, however, enhanced. Using partially conserved axial current one can show that the 
form factor gp is of order MN/niq ~ 100, making the contribution to the amplitude of order lO^"*. In 



Section 5.4 we review the status of experimental data and LQCD calculations showing this enhancement. 
The effect of gp on the neutron beta-decay rate has been worked out in Ref. |28j, and it should be 
included when the experiments reach that level of precision. 
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3.1 Differential decay distribution 

Including the effect of recoil corrections, radiative corrections, and BSM couplings, the 
differential decay rate for polarized neutrons reads [Ml ESI ESI EZ] 

Hr (f<W\2 \y \2 , . 



(7) 



where Pe and p^ denote the electron and neutrino three-momenta, while cr„ denotes the 
neutron polarization. The bulk of the electron spectrum is described by 

w{E,) = p,Ee{Eo-E,y F{Z = 1,E,) (^1 + + ^^i^H^e)) (8) 

where = A - (A^ - m^)/(2M„) (with A = M„ - Mp) is the electron endpoint 
energy, mg is the electron mass, and F{Z, Eg) is the Fermi function that captures the 
Coulomb radiative corrections {Z denotes the charge of the daughter nucleus, which 
coincides with the proton in this case). The function 5a\Ee) [Ml EZ] captures model- 
independent ("outer") radiative corrections, while the coupling Cy is sensitive to the 
short-distance ("inner") radiative correction [381 ES]. The differential decay distribution 
function D{Ee, Pe,Pu, Cn) is given by 



D(Ee,Pe,P.,^n) = 1 + Cq + C, + ^6 + a(i?e) + ^(E^) ^ 

+ mf-^ + c,...(Be) (5|^)' + c„.,(£.)5|^^ 

+ C,.s,{E.)^^, (9) 

where b is an effective Fierz interference term and a{Ee), A{Ee), B{Ee) and Caa,aA,aB{Ee) 
are effective energy-dependent correlation coefficients, whose full expressions [36l EZl 
[39] we report in Appendix |B| where one can also find the coefficients Co,i generated 
by recoil corrections]^ In absence of radiative corrections, recoil corrections and BSM 
contributions, the effective correlation coefficients d{Ee), A{Ee) and B{Ee) reduce to the 
following well-known leading-order expressions 

T/^^ 2A(1-A) 2A(1 + A) , , 

^(^J-ITsv- ^(^=>-tW' '1°) 

with the rest of coefficients (co,i, b, C( aa,aA,aB){Ee)) Vanishing in this limit. 
The impact of new-physics contributions can be summarized as follows: 

• The effect of e^/ij was already evident from the effective Lagrangian of Eq. |3j they 
induce (i) an overall correction proportional to (1 + 2eL + 2e/j), and (ii) the shift 
A — !■ A = A(l — 26^). As a consequence of this second effect, working to linear order 



See also Ref. [ID] for a discussion of recoil corrections to the proton asymmetry. 
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in new-physics contributions, the measurements of different correlation coefficients 
by themselves cannot disentangle A and eR, they simply provide independent mea- 
sures of A. In order to probe from correlation measurements, one needs to 
independently know Qa/ov from LQCD calculations. 

The scalar and tensor interactions es^r appear to linear order only through the 
Fierz interference term b and the analogue term in the neutrino-asymmetry 
parameter {bi^ is the part of B{Ee) proportional to nie/ E^,, see Appendix [b| for a 
precise definition) 



lBSM 



lBSM 



1 + 3A2 
2 

1 + 3A2 



gs es - 12A qt er 



OMgs es - 5.22 qtCt 



(11a) 



gs es\-4:gT (1 + 2A) 



0.44 gs es - 4.85 gx er • (Hb) 



To the order we are working, in the above expressions we can use either A or A. 

Experimentally, one can probe the new-physics contributions in A, 6^^^, and b^^^ 
through (i) measurements of the electron spectrum, aimed to isolate the term b in 
Eq. (|9|; or (ii) correlation measurements, aimed to isolate a{Ee), A{Ee), and B{Ee) 
in Eq. ([9]). Correlation measurements involve the construction of asymmetry ratios [35]. 
For example, in order to isolate A{Ee) one constructs the ratio A^xpiEe) = {N+{Ee) — 
N_{Ee)) / {N+{Ee) + N_{Ee)) where N±{Ee) are the spectra corresponding to events 
with cr„ ■ pe > and cr„ ■ pe < 0. Similarly, in order to isolate B{Ee) one can use the 
simple ratio 5exp(^e) = (Q++(^e) - Q—{Ee))/{Q++iEe) + Q—iEe)), where Q++(^e) 

and Q (i?e) are the spectra of events with cr„ ■ pe > 0, cr„ ■ pp > and cr„ ■ pe < 0, 

Cn ■ Pp < 0, respectively. One can immediately see that through the total spectra in the 
denominator, both y4exp(-Ee) and 5exp(-E'e) are sensitive to the Fierz interference term 
b, so that asymmetry measurements involving simple ratios as described above really 
measure 

yiEe) 



Y{E, 



(12) 



1 + bnif,/ Ef. 

where Y G {A, 5, a, ...}. Moreover, each individual experiment applies optimization cuts 



in Ef,, thus mesuring a specific weighted average of Eq. (12). 



The above observation has important consequences for the phenomenology of neutron 
decay: (i) The rrie/Ee component of Bexp{Ee) is sensitive not to b^^^ but rather to the 
combination (1 + 3A2)/(2A(1 + A))6f ^ - b^^^ ^ - b^^^. Besides 5exp(^e), it 
might be possible to construct a set of observables that disentangle the contribution of 
^BSM ^BSM ^Yiis case the BSM sensitivity of b^^^ alone is of interest. In 

our phenomenological analysis we will study both cases (constraints from b^^^ — b^^^ 
and b^^^). (ii) More generally, correlation coefficients measurements traditionally used 
to determine A = gA/gv within the SM {ei/R = 0,b = b^ = 0), provide information 
on three independent parameters in our EFT setup: A = A(l — 2e/j), b^^^, and b^^^^^ 



^ In other words, if es.T are larger than the experimental errors, one has to observe an unexpected 
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A fit to the current data |l2l iSl iU gS] (with precision 5A/A ~ 0.005, 5a/ a ~ 0.05, 
5B/B ~ 0.005) yields -0.3 < 6^sm j,bsm ^ g 5 95% q l g], which, as we will 

see, is not competitive with other bounds. It will be interesting, however, to explore the 
implications of future experimental improvements in the combined extraction of A, 5^^^ 
and 6^^^ from a. A, and B measurements, along the lines described in Refs. [H H6] . 

The main conclusion from the above discussion is that measurements of the differ- 
ential neutron-decay distribution are mostly sensitive to new physics through 5^^^ and 
jjBSM^ which depend on the scalar and tensor couplings, and e^, to linear order. There- 
fore, apart from the next section, which we include for completeness, in the rest of this 
paper we restrict our discussion on these exotic scalar and tensor interactions, comparing 
the physics reach of neutron decay to other low-energy and collider probes. 



3.2 Total decay rate and determination of 



ud 



For completeness, we discuss here the BSM corrections to the neutron decay rate and 
the extraction of Vud from neutron decay. Expressing in terms of the Fermi constant 
determined in muon decay (this involves non-standard contributions to the purely 
leptonic charged-current interaction encoded in the coefficient [21]) and performing 
the phase-space integrations, the total decay rate reads 



Gl\Vud\^m 



27r3 



I + 3P ■/■ (1 + A 



RCj 



1 + 2eL - 2vl + 2eR + h 



BSM 



JiOro) 



■ (13) 



In the above expression, the corrections from BSM physics are encoded in A and the terms 
in square brackets. Arc = 3.90(8) x 10~^ is the SM electroweak radiative correction [38j, 
and the phase-space integrals are defined by 



^ (xo 



X 



)Vx2 



1 dx 



/ = /o(xo)(l + Ay) 



(14) 



where xq = Eo/rrie and Af encodes Coulomb and recoil corrections that are numerically 
quite important, /o(a;o) = 1.629, / = 1.6887, Ji(xo)//o(a;o) = 0.652 (See Ref. [38] for 
details). In order to extract Vud from neutron decays one needs (see Eq. [l3| experimental 
input on the neutron lifetime l/F [HI US] and A, which is usually extracted from beta- 
asymmetry Ae:^p{Ee) measurements [^ H3] (after accounting for recoil and radiative 
corrections). Taking into account Eq. (12), the usual method for extracting A actually 



determines A (l + c6^^^), where c is a certain 0(1) number that depends on the specific 
experimental analysis. In summary what we really extract from neutron beta decay is 

energy dependence of the form m/E in the measurements of the correlation coefhcients (in addition 
to the various expected energy dependences due to sub-leading standard effects that are detailed in 
Appendix [b]) . Thus, for a certain energy, a determination of A from a{A) would be actually extract- 
ing the quantity X (l + n^j^j^^b^^^m/ , whereas in a i3-based determination of A, we would have 
A (1 + nB(6^S^ - b^s^)m/E), where = i^-^'Ki+sx') ^ _q 28, ua - « -0-25 and 

(l+A)(l+3A^) ^ 2 . 



(1-A)(1+3A) 
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not Vud but the combination 





= \Vud\^ 


l + 2eL 


-2VL 


+ 26^ + 6^^^ 
















l + 2eL 


-2VL 


+ 2e^ + 6^s^ 



/o(a;o) 1 + 3A2 



(15) 



4 Low-energy phenomenology of scalar and tensor 
interactions 

4.1 Other probes of scalar and tensor interactions 

In order to assess the discovery potential of experiments planning to measure h and B 
at the level of 10~^ and 10~^, it is crucial to identify existing constraints on new scalar 
and tensor operators. As we discuss below in some detail, the most stringent constraint 
on the scalar coupling arises from O"*" — O"*" nuclear beta decays. On the other hand, 
the most stringent bound on the tensor effective coupling arises from the Dalitz-plot 
study of the radiative pion decay tt — t- ez/7. For completeness, we will also briefly review 
(i) constraints on es^r from other nuclear beta-decay observables, showing that they are 
not competitive at the moment; and (ii) constraints on es^p^r arising from the helicity- 
suppressed tt — )• ez/ decay. As we will show, the latter provides potentially the strongest 
constraints on es^r-, once the flavor structure of the underlying theory is known. This 
provides very stringent constraints on model building. 

4.1.1 0^ — i- 0"*" transitions and scalar interactions 

At leading order within the SM and new physics, the differential decay rate for an 
unpolarized nucleus is [32] 

'^^'"^'^ - ^ (1 + 26, + 2en) PeE.{h - E^fF{-Z^ E,) 



dEedVLedVty (2ti] 
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where Eq = Mp — Mp, is the electron endpoint energy expressed in terms of the masses 
of parent and daughter nuclei, F{—Z, Ef.) is the Fermi function, Z is the atomic number 
of the daughter nucleus (the minus sign applies to emitters for which the most precise 
measurements exist). For O"*" — )■ 0"^ transitions the coefficients a, 6 are 

ao+ = 1 (17a) 
60+ = -2igses 7 = - a^Z^ , (17b) 

and the total rate is given by 



_ Gl\Vud?rnl / (0+^0+) 

o+-s>o+ — /o+-H>o+ y- + 



l + 2ez.-2{}i + 2e^ + 6o+-^'^^°^ 



;i8) 
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where xq = Eo/nie. In this last expression, the SM sub-effects have been included 
through aII'q"^'^ ^ and also inside /o+_s>o+, that up to Coulomb, nuclear distortion and 
recoil effects, is /o+_s>o+ = loi^o), similarly to what happens in the neutron-decay case. 
The various radiative corrections (including Ajj'^^"^'^ ■*) are discussed in detail in Refs. [IHl 
150] . Comparing the values of Vud as extracted from neutron and nuclear decays, we find 



see Eq. ( 15 ) and the preceding discussion) 



- 1 + - 6f ^ ( ^44 - ) , (19) 



which in principle provides another handle on scalar and tensor interactions. 

Let us now come to the point of greatest interest for this paper's discussion. From a 
comparison of precisely known half-lives corrected by phase-space factors /o+_!.o+, Hardy 
and Towner |50j found 6o+ = —0.0022(26), which translates into the following bound on 
the product of nucleon scalar form factor and short-distance scalar coupling: 

- 1.0 X 10"^ < ^5 65 < 3.2 X 10^3 (90% C.L.) . (20) 

This is the most stringent bound on scalar interactions from low-energy probes. 



4.1.2 Radiative pion decay and the tensor interaction 

An analysis of the Dalitz plot of the radiative pion decay vr"*" — )• e'^u^'j is sensitive to 
the same tensor operator that can be probed in beta decays. The experimental results 
from the PIBETA collaboration [51] put constraints on the product ex x fx of the short- 
distance coupling and the hadronic form factor /r defined by [52] 

{l{^,p)\u(T^ul5d\7r^) = "2 /t (Pi^^u-Pue^) , (21) 

where and are the photon four-momentum and polarization vector, respectively. 
The analysis of Ref . [S2] , based on a large- A^c- inspired resonance-saturation model pro- 
vides /t = 0.24(4) at the renormalization scale /z = 1 GeV, with parametric uncertainty 
induced by the uncertainty in the quark condensate. The 90%-C.L. experimental con- 
straint]^ —2.0 X 10"'' < er X /t < 2.6 x 10~^, when combined with the above estimate 
for fx run to 2 GeV implies 

- 1.1 X 10-3 < er < 1.36 x 10"^ (90% C.L.) . (22) 

Again, this is the most stringent constraint on the tensor coupling from low-energy 
experiments. The next best constraints, which we report in the next section, arise from 
measurements of nuclear beta decays. 

^Note that there is a factor of 2 difference in the normalization of the tensor couphng er compared 
to what was used in Refs. [25l I51j . 
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4.1.3 Bounds on scalar and tensor structures from other nuclear beta decays 

Bounds on scalar and tensor interactions can be obtained from a number of observables 
in nuclear beta decays, other than O"*" — ?• O"*" transitions. Although these bounds are 
currently not competitive, we summarize them here for completeness. 

The leading sensitivity to scalar and tensor operators appears through the Fierz 
interference term b, which in the limit of pure Gamow- Teller transitions is proportional to 
the tensor coupling (6gt = —{^'JOt^t) / while in pure Fermi transitions is proportional 
to the scalar coupling [bp = 2'-fgs€s)- Significant constraints on b arise from electron- 
polarization observables [32] as well as in measurements of A and a in both Fermi and 
Gamow- Teller transitions. Here is a summary of current bounds on es,T'- 

• The most stringent constraint from the beta asymmetry in pure Gamow- Teller 
transitions (Aqt) arises from ^"Co measurements and implies [53] 

-2.9 X 10"^ < grer < 1.5 x 10"^ (90% C.L.) . (23) 

Similar bounds can be obtained from measurements of Aqt in ^^^In decay |54j : 
-2.2x10-2 < g^er < 1.3 x 10"^ (90 % C.L.). 

• Measurements of the ratio Pf/Pqt of longitudinal polarization in the positron 
emitted in pure Fermi and Gamow- Teller transitions [55l 156] imply 

4 

- 0.76 X 10-2 < gses + -^grer < 1-0 x lO'^ (90% C.L.) . (24) 

A 

• Preliminary results have been reported on the measurement of the longitudinal 
polarization of positrons emitted by polarized ^°^In nuclei [57] . The corresponding 
90 % C.L. sensitivity to tensor interactions, Igx^r] < 3.1 x 10"^, is quite promising 
although not yet competitive with the radiative pion decay. 

• Finally, the beta-neutrino correlation a has been measured in a number of nu- 
clear transitions [581 EH EO] [61]. The resulting constraints on scalar and tensor 
interactions are nicely summarized in Fig. 7 of Ref. [58]. In terms of the coupling 
constants used here, the 90 % C.L. combined bound on the tensor interaction reads 
Ifi'T^rl < 5 X 10"^, again not competitive with the radiative pion decay. 

We observe that in order to improve on the existing bound on from vr — )■ ewy, future 
measurements sensitive to 6gt should aim at sensitivities of 5&gt ^ 6.3 x gT>^ 10^^. For 
example, a 10~^ measurement of 6gt would probe gr^T at the 2 x 10~^-level, providing 
a very competitive bound. 

4.1.4 Constraints on es^p.T from vr — ez/ 

The ratio Rt, = T^tc — )■ eh'['y])/T{n — /iz/[7]) probes more than just the effective low- 
energy pseudoscalar coupling ep defined earlier as the coefficient of the operator e(l — 
75)z/e ■ U'j^d. In fact, since (i) i?^ is deffned as the ratio of electron-to-muon decay and 
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(ii) the neutrino flavor in both the decays is not observed, this observable is sensitive to 
the whole set of parameters €p^ defined by 



eff 



D 



V2 



T A C 

Vud e; 



,(1 - 75)z/^ ■ u-f^d 



(25) 



where a G {e, fi} refers to the flavor of the charged lepton and /3 G {e, yU, r} refers to 
the neutrino flavor. One generically expects SM extensions to generate non-diagonal 
components in e^^j,^ In the new notation the previously defined pseudoscalar, scalar, 
and tensor couplings reads ep^s,T = ^pst- It is important to note here that only ep and 



can interfere with the SM amplitudes, while the remaining ep^ contribute incoherently 
to both the numerator and denominator in -RttQ In summary, allowing for non-standard 
interactions and factoring out the SM prediction for i^^r, one can writ^ 



-Bo ^ee 



+ 



Bo ef^ 



Bo, 



rriu P 



Bp ^IJ.e 



Bp ^f^T 
rriu P 



(26) 



In the above equation the factors of BQ/nie^^ep represent the ratio of new-physics am- 
plitude over SM amplitude. The latter is proportional to the charged-lepton mass 
due to angular-momentum conservation arguments, while the former is proportional 
to (0|'U75(i|7r), characterized by the scale- and scheme-dependent parameteip 



Ml 



(27) 



Since B^^{n = 1 GeV) = 1.85 GeV and consequently Bo/rUe = 3.6 x 10^, i?^ has 
enhanced sensitivity to ep^, and one needs to keep quadratic terms in these new physics 



coefficients 



Inspection of Eq. (26) reveals that if the new-physics couplings respect e^/mg 



ep°/m^, then Rt^/E^ = 1, and there are no constraints on these couplings. On the 
other hand, if the effective couplings ep'^ are all of similar size, one can neglect the 
entire denominator in Eq. (26), a it is suppressed with respect to the numerator by 



powers of mg/m^. We will assume to be in this second scenario. In this case the con- 
to live in a spherical shell of radius 



eri 



er 

P 5 '^P 



straint in Eq. (26) forces the couplings ep 

mjBo^RT^/R^ ^ 2.75 x 10"^ centered at e'jS = mjB^ ^ 2.75 x 10-^ e 



0. 



^ While in our setup the incoherent contribution arises from "wrong-flavor" neutrinos, in general it 
could have a different nature. For example, the incoherent contribution to discussed in Refs. [551 [20] 
is due to a right-handed light neutrino. 

^Here we are neglecting the overall effect of Ui/fl, not enhanced by helicity arguments. 

^Note that the scale and scheme dependence of -Bo(m) is compensated in physical quantities by the 
scale and scheme dependence of the Wilson coefficients ep^ . 

^'^ This feature is specific to purely leptonic decays of pseudoscalar mesons. In beta decays one never 
encounters relative enhancement factors such as B/rrie, because ep is always multiplied by nucleon 
velocity factors and the SM amplitude does not suffer anomalous suppression (as the helicity argument 
implies in the case of tt — > ev) . 
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Figure 1: The allowed region in the two-dimensional plane 
given by an annulus of thickness 1.38 x 10 
90 % C.L. bound on e!? is -1.4 x 10 



ep-ep determined by R^, is 
In the absence of information on eff, the 



< < 5.5 X 10" 



The thickness of the shell is numerically 1.38 x 10 ^ and is determined by the current com- 
bined uncertainty in R"^^ [S21 E3] and Rf^ [BllEl]: R^^^/Rf^ = 0.996(5) (90% C.L.). 
This is illustrated in Fig. [T| where we plot the allowed region in the two-dimensional 
plane given by ef? and a generic "wrong- flavor" coupling denoted by ef? . Note that the 
allowed region is given by the thickness of the curve in the figure, thus enforcing a strong 
correlation between ep and ep . Since e'^^ are essentially unconstrained by other mea- 
surements and can be of order 10~^, we can marginalize over either one of the couplings 
to obtain a bound on the other. The resulting 90%-C.L. bounds are 



-1.4x10"^ < e"^ < 5.5x10' 



or 



-2.75x10"* < < 2.75x10" 



(« 7^ e) , (28) 



in qualitative agreement with the findings of Refs. [25l 120] . 

As originally discussed in Refs. [HI [25l |26], the pseudoscalar coupling cp 
radiatively generated starting from nonzero es,T- Hence, the stringent constraint in 



elf can be 



Eq. (28 ) puts constraints on the same e^^^ that can be probed in beta decays. The physics 
of this effect is very simple: once the scalar, pseudoscalar, and tensor operators are 
generated by some non-standard physics at the matching scale A, electroweak radiative 
corrections induce mixing among these three operators. So even if one engineers a small 
pseudoscalar contribution ep(A) at the matching scale, known SM physics generates a 
nonzero ep(/i) at some lower energy scale /i via loop diagrams. The general form of the 
constraint can be worked out by using the three-operator mixing results from Ref. p^^^"^ 



'^Thc authors of Ref. [2_6 focused only on the phenomenology of scalar-to-pseudoscalar mixing. 
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The leading-order result is 



ISP 
iTP 



(29a) 



ep'^(A) ( 1 + 7pp log - j + e5^(A) 7sp log - + e^'^(A) 7tp log - 

113 ttl n , , s 

1.3 X 10^2 (29b) 

(29c) 
(29d) 



3^2 

4 vr 72 vr 

15 a^ . 
^ ^ 6.7 X 10"^ 



72 TT 
9 OL2 

^2 V 



15 a\ n 
y— ^ -7.3 X 10^2 



TT 



where a\ = a/ cos^6'vi/ and ^2 = «/ sin^ ^v(/ are the U{1) and SU{2) weak couplings, 
expressed in terms of the fine-structure constant and the weak mixing angle. Setting 
ep(A) = and neglecting the small 0{a/Ti) fractional difference between es^ri^) and 

{) C.L. constraint on the es-er plane reads 

5.5 X 10"^ 



the observable es,T{lj) at the low scale, the 
-1.4 X 10"^ 



log(A//i) 



< ISP £5 + iTP Ct < 



log(A//i) 



(30) 



Even assuming log(A//i) ~ 10 (e.g. A ~ lOTeV and /i ~ IGeV), using the numerical 
values of 7s'p,tp, one can verify that the individual constraints are at the level of \es\ ^ 
8 X 10~^ and le-rl < 10~^, implying that this constraint on e-r is roughly equivalent to 
the one arising from vr — > ez/7. Of course, these bounds become logarithmically more 
stringent as the new-physics scale A grows. 



to 
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Figure 2: 90% C.L. allowed regions in the e^-e^ plane implied by (i) the existing bound 
on 60+ (green horizontal band); (ii) projected 10~^-level limits on b (red band), — b 
(blue band, left panel), and b^ (blue band, right panel). The hadronic form factors are 
taken to he gs = Qt = ^ in the ideal scenario of no uncertainty. The impact of hadronic 
uncertainties is discussed in Section |6l 
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Figure 3: 90% C.L. allowed regions in the e^-er plane implied by (i) the existing bound 
on 60+ (green horizontal band); (ii) projected 10~^- level limits on b (red band), by — h 
(blue band, left panel), and by (blue band, right panel). The hadronic form factors are 
taken to be (75 = = 1 in the ideal scenario of no uncertainty. The impact of hadronic 
uncertainties is discussed in Section |6l 



4.2 The impact of future h and B neutron measurements 

The discussion in the preceding subsection has shown that currently the most stringent 
low-energy constraints on novel scalar and tensor interactions arise, respectively, from the 



Fierz interference term in O"*" — )■ O"*" nuclear beta decays (Eq. (20 )) and from the radiative 



pion decay tt — > ez/7 (Eq. (22)). It is important to realize that the allowed interval 
derived from Eq. (20) depends on the nucleon form factor gs (as do all the constraints 
arising from neutron and nuclear beta decays). For a given experimental accuracy, the 
constraint on the short-distance couplings es,T becomes stronger as Sgs,T/gs,T — ^ 0. In 
this section, we will first explore the maximal constraining power of nuclear and neutron 
measurements in the ideal scenario of no uncertainty on gs,T, and for illustrative purposes 
we assume the central values gs = gr = ^- We will quantify the implications of finite 
uncertainties on gs^^ on the es,T constraints in Section [6j 

With the above assumptions on gs,T, the currently allowed region (at 90% C.L.) on 
the es-^T plane is given by the green horizontal band in Figs. |2]and[3] The vertical {es) 
boundaries of this region are determined by the constraint from 69+ 1 while essentially the 
entire horizontal (e^) range on the scale of these plots is allowed by the tt — )■ ez/7 limit 



(see Eq. ([22])). 

In this ideal scenario of no uncertainty on gs,Ti we can quantify the impact of future 
neutron measurements by plotting the 90% C.L. allowed region in the e^-e^ plane implied 
by projected limits on 6, by — b, and by. The neutron constraints are derived using 



Eqs. (11) and in generating the plots we use the central value A = 1.269. In Fig. [2] we 
focus on the case in which the experimental sensitivity on b, by — b, and by is at the 
10~^ level. In the left panel we show the constraints from the existing 60+ limit (green 
horizontal band) and 10~^-level limits on b and by — b, (red and blue bands, respectively). 
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In the right panel we replace the 10~^-level limit on — b with the 10~^ limit on b^, 
which in principle can be isolated experimentally [JT]. In Fig. [sjwe plot the constraints 
resulting from projected limits on b, b^ — 6, and by at the 10~^ level. The intersection 
of the various bands in Figs. [2] and [3] denotes the combined allowed region in the e^-e^ 
plane that would result after future neutron measurements. Two important remarks are 
in order here: 

• For a given experimental sensitivity, the combination by — b gives weaker constraints 



on es^T than b or by. This is easily understood: by taking the difference of Eqs. ( 11 ) 
one sees that by — b (x \ — 1, which for A ~ 1.27 provides a suppression factor. 

There is an almost exact "degeneracy" in the constraints from b and by, again 



controlled by the form of Eqs. (11) and the numerical value of A. For the purposes 
of constraining es^x-, an upper limit on b is essentially equivalent to an upper limit 
on by. This provides strong motivation to pursue experimental determinations of 
both by — b and by via neutrino asymmetry {B) measurements. From the theoretical 
point of view, we can use either b oi by, and in subsequent sections we will use b 
for illustrative purposes. 

Fig. [2] clearly illustrates that with experimental sensitivity in neutron decay at the 



'0+ 



10~ level, the most stringent constraint arises from a combination of 6o+ and b or 6| 
and by. The complementarity of these measurements would lead to a significant (four- 
fold) improvement in the bound on ey, compared to Eq. (22). The impact of lO"'^ 



measurements of 6, by, and by — b'm. neutron decay is even more dramatic (Fig. [3]), as in 
that case the constraint from 6o+ would become irrelevant and the combination of b and 
by — b or b and by would imply an improvement of one order of magnitude in the bound 
on e-r and a factor of two in e^. 

In Section |6] we will revisit the impact of proposed neutron measurements on es,T in 
light of nonzero uncertainties in the hadronic matrix elements gs,T- 



5 Lattice calculation of matrix elements 

To connect the measurements of b and by in neutron decays to new physics at the TeV 
scale requires precision measurements of the matrix elements of isovector bilinear quark 
operators between an initial neutron and final proton state, in particular of the scalar 
and tensor operators. Lattice QCD is a path-integral formulation of QCD on a discrete, 
four-dimensional Euclidean spacetime, and numerical simulations of it provide the best 
nonperturbative method for evaluating these matrix elements. It has been successfully 
employed to calculate hadron masses and their decay properties, such as matrix elements, 
with control over statistical and all systematic errors, in many cases at higher precision 
than can be measured experimentally [661 EZj • 

To obtain continuum results, estimates from LQCD obtained at a number of values 
of lattice spacing a and spacetime volume xT are extrapolated to a — )■ and L — )■ oo 
to eliminate the artifacts introduced by formulating QCD in a finite discretized box. An- 
other source of systematic uncertainty is introduced when estimates obtained at multiple 
values of u and d quark masses heavier than in nature are extrapolated to the physical 
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point. One typically uses cliiral perturbation theory to carry out this extrapolation, 
with low-energy constants determined by over-constraining the fits using experimental 
and lattice data [HH]- Current state-of-the-art simulations are beginning to provide re- 
sults at physical light-quark masses obviating the need for a chiral extrapolation. Recent 
calculations by the BMW collaboration [69], [70] at multiple lattice spacings, volumes and 
pion masses as light as 120 MeV provide an excellent demonstration of how hadronic 
properties can be extracted with fully understood and controlled systematics. 

In this section we review current LQCD calculations of the nucleon isovector matrix 
elements in order to highlight what needs to be done to obtain the precision required to 
probe new physics at the TeV scale in neutron-decay experiments. We also present our 
current best estimates of gs and qt, which are used in the phenomenological analysis 
presented in Sec. |6] 

5.1 Lattice methodology 

A lattice calculation proceeds in two steps: First, a Monte-Carlo sampling of the QCD 
vacuum, called an "ensemble of gauge-field configurations", is generated using an ap- 
propriate discretization of the gauge and fermion actions. The particular choices of the 
actions have important implications for the computational cost of the calculation, for 
the size of the discretization errors and for which symmetries are violated at finite lattice 
spacing. We will review the existing calculations, summarized in Table [T| with two light 
flavors (2-fiavor) and two light and one strange flavor (2-|-l-fiavor) as these are close 
approximations to the real world. 

The second step is to calculate expectation values on these ensembles of gauge con- 
figurations and from these extract estimates of the desired observables. For hadronic 
observables, the fermion action used at this stage may differ from the one used in 
producing the gauge configurations, in which case it is called a "mixed-action" cal- 
culation. Further details on the domain-wall fermion (DWF) formulation are given in 
Refs. [m [721 ESI EH ESj ; clover fermions in Ref. [TB]; twisted- mass fermions in Ref. [77] : 
and improved staggered fermions in Refs. [TSl [721 [SO] . 

Calculation of the isovector nuclear matrix elements requires two separate optimiza- 
tions in addition to the choice of the actions. The first is to tune the size of smearing 
applied to the local interpolating operator with the correct quantum numbers of the 
nucleon 

X'^ix) = eatc[^fix)C^,^p',ix)]riix), (31) 

where a, b, c are color indices, C is the charge-conjugation matrix, and ipi and il'2 
are u or d quarks; for example, to create a proton, we want ipi = u and %l)2 = d. 
This local operator, unfortunately, couples to the nucleon and all its excited states with 
the same quantum numbers. To improve the overlap with the desired ground state, 
the quark fields in this operator may be "smeared" around the point x. The goal of 
this smearing is to approximate the ground-state nucleon wavefunction. We adopt the 
commonly used application of the three-dimensional gauge-invariant Laplacian to smear 
around the source point x and tune the smearing size to improve the overlap with the 
ground-state nucleon in the two- and three-point correlation functions. The two-point 
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function, projected to a definite momentum at either the source or sink time by making 
a three-dimensional Fourier transformation, is given by 

rS(t;p) = {x^{t,p){x^np)) 

- E(0|XA(^,p)|r^>(n|(x^)^(0)|0>^^e-^"(P)*, (32) 

where the indices A and B indicate the choice of operator smearing. The nucleon states 
are normahzed as (0|(x;A)^b, s) = XaUs{p) where Xa is the overlap of the operator with 
the state, and the spinors satisfy Ms(p)its(p) = -E'(p)7* — ■ p + m. In the limit 
of large time separation t, the correlator is dominated by the ground-state nucleon, and 
the above form simplifies to 

rS«;p) = ^?^|ipx.(p)x,(p)e--ip)'. (33) 

To calculate the nucleon matrix elements, we also need to construct nucleon three- 
point functions with insertion operators Oy{x) = Z^O^ — ZYu{x)Td{x), where is the 
bare operator, F represents one of the sixteen Dirac matrices and Zr is the associated 
renormalization constant of the operator. The three-point functions take the form 



tf- p„ p;) = Zr {^x% itf, PfWmi iU: P^)) ■ (34) 

By inserting a complete set of states {n, n'} between the operators, this three-point 
function can be written as 

i^, (t.,M/,P.,P/) AE'^{pf)EM 

X ^T^pUn'{Pf,s')fj{Nn>{pf,s') \O^Nn{pi,s))Un{Pi,s)a, 
s,s' 

(35) 

where T is an appropriate projection on the baryon spinors. At sufficiently large source- 
sink separation (tgep = tf — ti), the signal due to excited states dies out exponentially, and 
the sum over states reduces to just the ground states n = n' = 0. The operator overlap 
factors Xa,b and the exponential time dependence can be canceled out by constructing 
a ratio of three- and two-point functions, which for the simple case of pj = p/ = is 



rf /{U,t,tf;p, = 0,pf ^ 0) 
rf/{ti,tf;p^O) 



= ^(2),T,. . :^ ■ (36) 



In practice, choosing a sufficiently large source-sink separation tgep to make the excited- 
state contamination negligible is challenging because the statistical signal in both the 
two- and three-point functions involving nucleons degrades exponentially with tgcp- Thus, 
the second optimization required is over tgep- In ongoing LQCD calculations we are 
exploring multiple values of igep and will explicitly include excited states in our analysis 
to understand and reduce this systematic error. 
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Collaboration 


Action 




A-U (McV) 


L (fm) 




n a (fm) 


Calculated 


QCDSF|8ll 


clover 


2 


595-1000 


1.0-2.0 


4.6 


0.07-0.116 


9 A 


QCDSFina 


clover 


2 


170-270 


2.1-3.0 


2.6 


0.08-0.116 


9 A, 9T 


CLS[83] 


clover 


2 


290-575 


1.7-3.4 


4.2 


{0.05,0.07,0.08} 


9A 


ETMC[84] 


twisted Wilson 


2 


260-470 


{2.1, 2.8} 


3.3 


{0.056,0.070,0.089} 


9A 


RBCSg 


DWF 


2 


490-695 


1.9 


4.75 


0.117 


9A, 9p, 9t, 9v 


RBC/UKQCD[8l[87j 


DWF 


2+1 


330-670 


{1.8,2.7} 


3.8 


0.114 


9 A, 9t 


LHPCIHSIIMIEO] 


DWF on staggered 2+1 


290-870 


{2.5,2.7} 


3.68 


0.1224 


9A, 9*p, 9T 


QCDSF[91] 


clover 


2+1 


350-480 


1.87 


3.37 


0.078 


9A 


HSC,92] 


anisotropic clover 2+1 


450-840 


2.0 


4.57 


0.125 {at = 0.036) 


9A 



Table 1: A summary of recent LQCD calculations of qa, g*p and gx by different collab- 
orations using two and three flavors of dynamical quarks and 0(a)-improved actions. 



For brevity, we use g*p for the induced-pseudoscalar charge discussed in Sec. 5.4 and 
(M,r-^)min for the minimum value of Mj^L used in that set of calculations. 



5.2 Issues in extracting the matrix elements 

The matrix elements of most interest to us are those of the scalar and tensor bilinear 
operators, ud and ua^„d] however, we are calculating all five Lorentz structures as the 
additional cost is negligible. There is independent interest in high-precision measure- 
ments of gAi and it provides a cross-check of the lattice systematics. The three-point 
correlation functions of the vector operator will be used to construct ratios of matrix 
elements and renormalization constants to reduce systematic errors. In this section we 
summarize issues relevant to the LQCD calculations of these matrix elements. 

The desired matrix elements of isovector bilinear operators Oy{x) = Zru{x)rd{x) 
have a number of simplifying features and allow us to make certain approximations: 

• There are no disconnected Feynman diagrams contributing to the three-point func- 
tions. These typically arise when quark fields in composite operators can be con- 
tracted between themselves. 

• There are no lower-dimensional operators with which isovector bilinear operators 
mix, so there are no power-law divergences. Only multiplicative renormalization 
factors Zr need to be calculated. 

• Current lattice simulations are done with degenerate u, d quarks, at zero momen- 
tum transfer, and do not include electromagnetic effects. The momentum transfer 
in neutron decay, = 1.7 MeV^ is sufficiently small that the matrix elements can 
be calculated at = 0. Also, the isospin-breaking and electromagnetic contribu- 
tions are expected to be smaller than the statistical errors. 

• Protons and neutrons are both stable asymptotic states of strong interactions, so 
there are no other hadronic final states that complicate the calculations. 

The issues that need to be addressed to obtain precision results are the following: 

• The signal-to-noise ratio in both two- and three-point correlators decreases rapidly 



with the time separation tgep between the source and the sink in Eq. 36 It is 



therefore, necessary to improve the signal by increasing the overlap of the operators 
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used as sources/sinks with the nucleoli ground state. As discussed in Sec. 5.1, our 



current approach is to (i) smear the quark fields in the interpolating operator given 



in Eq. 31 and tune the smearing size, and (ii) explicitly include excited states in 



the analysis. 

A careful optimization of the Eucledian time interval tsep between the source and 
sink in the three-point functions has to be carried out for each lattice spacing a 
and light-quark mass. On the one hand, this interval should be as large as possible 
to isolate the nucleon ground state on either side of the operator insertion, and 
on the other hand the statistical noise limits the time separation. While there is 
no a priori minimum value of tsep as it depends on how well the source and sink 



operators are tuned, in Sec. |5.3| we show that current data suggest that asymptotic 
estimates are obtained with tgep > 1.2 fm for the operators used. Our focus will 
be on improving the operators and investigating 2-3 values of tsep to reduce and 
quantify this systematic error. 

• One needs to demonstrate that the lattices are large enough that finite-size effects 
are under control, especially for proposed calculations with pions masses below 
350 MeV. When the spatial volume used is too small, finite- volume effects arise due 
to the coarseness of the available lattice momenta, squeezing of the wavefunction 
due to the interaction of a spatially extended particle with itself and contamination 
from partons wrapping around the lattice. Previous studies have shown, as a rule 
of thumb, that finite-size effects are smaller than statistical errors for M^L > 4. 
The detailed form of the finite- volume corrections is quantity-dependent. 

• Very high-statistics measurements, typically on a few thousand gauge configura- 
tions, will be needed to improve the signal in the two- and three-point correlation 
functions to overcome the rapid growth in noise with tggp- Our ongoing calcula- 
tions show that the statistics needed will be determined hy gs as it has the smallest 
signal-to-noise ratio. 

• The calculations need to be performed at a sufficient number of values of the light- 
quark mass to extrapolate results to the physical value mi = 0.037ms, and at 
sufficient number of values of the lattice spacing a to extrapolate to the continuum 
limit. 

• The renomalization constants Zr depend on the choice of both the gauge and 
fermion actions and have to be calculated for each ensemble of gauge configurations. 
In past calculations, Za typically varied between 0.75-0.9 for the lattice spacings 
that have been simulated. The scale-dependent Zs, Zp, and Zt (given in the MS 
scheme at 2 GeV) show larger variations and dependence on the lattice action. 
One-loop tadpole improved perturbation theory can underestimate corrections to 
|1 — Zr| by 50%. Nonperturbative methods, such as calculating in the RI-MOM 
scheme [021 EH ES], are preferred as they reduce this uncertainty to a few percent, 
and we will use them in our calculations. 

In the next four sub-sections we summarize the extent to which these issues are under 
control in current calculations of each of the matrix elements in order to highlight what 
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needs to be done to achieve the desired precision of 10-20%. The analyses of qa, Qp 
and qt are reviews of existing calculations, and the new estimate of gs we present is 
preliminary. 



5.3 Nucleon axial charge qa 

The axial charge of the nucleon qa = gA{<f' = 0), defined in Eq. [sj is a fundamental 
hadronic observable, well measured in neutron beta-decay experiments: qa = 1.2695(29) x 
9v [SS] , where the vector charge gy = 1 since Vud has been factored out in the Lagrangian 
given in Eq. ([s]). Since the axial charge is experimentally well known, it has long served 
as a benchmark quantity for LQCD calculations, particularly for estimating systematic 
errors in other nucleon matrix elements that are either poorly measured in experiments 
or completely unknown. 

Many groups worldwide have calculated qa using various gauge ensembles and fermion 
actions as summarized in Table [T] and shown in the first two panels of Fig. |4] for two- 
and three-flavor simulations, respectively. The results from each study, after a chiral 
extrapolation to the physical pion mass, are shown in the third panel of Fig. |4} 

The overall observations are: (i) The central values vary between 1.12 < qa < 1.26, 
and the errors are much larger than the experimental uncertainty. The deviations from 
the experimental value are large, considering that corrections due to strong interactions 
determine qa — ^- (ii) There is no significant difference between 2- and 2+1-fiavor 
estimates or dependence on the light-quark mass at these unphysically large oc m^. 
More high-precision calculations are needed to determine whether the chiral behavior 
changes at smaller quark masses and to gain control over the extrapolation to the physical 
M^. (iii) Within errors, the lattice data are consistent between the different groups 
(with different lattice actions), different lattice spacings and between 2- and 2 -|- 1-fiavor 



theories. Our understanding of systematic errors, discussed in Sec. |5.2[ are summarized 
next. 

Investigations of finite- volume effects have been carried out by the RBC/UKQCD 
collaboration [Sn]. They used domain- wall fermions at a fixed lattice spacing of 1/a = 
1.73(3) GeV (equivalently, a = 0.114(2) fm) on two lattice sizes L = 1.8 and 2.7 fm. 
They found that at fixed ~ 0.1 GeV^ there are significant finite-volume effects for 
L < 2.5 fm, and these lower the value of qa- They also analyzed qa as a function of 
M^L and found that the data scale in this variable; i.e. data from a given action and for 
a given number of flavors collapse onto a single curve. For small M^L, the value of qa is 
underestimated and to get within 1% of the infinite- volume result requires M^L > 6. 

The QCDSF collaboration [81] analyzed qa at four lattice spacings ranging from 
0.07 to 0.116 fm, and found no significant dependence on the lattice spacing. They, 
and the ETMC collaboration [84\, have also analyzed their data using finite- volume 
corrections suggested by heavy-baryon chiral perturbation theory (HBxPT) with small- 
scale expansion. They find that correcting their data for finite-volume effects at each 
lattice spacing improves their extrapolation to the physical pion mass. On the other 
hand, the RBC collaboration [85] finds that such corrections do not account for their 
data either qualitatively or quantitatively. An understanding of finite-volume effects, 
therefore, needs more work. 



24 



A source of potentially large systematic error is excited-state contamination when the 
source-sink separation (tgep = if — ti) is insufficient. The 2008 RBC 2- flavor study [85] 
used t ^ l.Q and 1.2 fm to check whether there is any significant dependence of qa on 
tsep- In this case, the values of qa were consistent within statistical error, and the central 
values increase by less than 5% between tgep = 1-0 and 1.2 fm. Similarly, LHPC [90] 
observe a tiny shift of the central value when changing tgep from 1.1 to 1.2 fm. The 
recent work by the CLS collaboration [83] investigated smaller source-sink separations, 
^scp ~ 0.56, 0.70 and 1.05 fm, on their a ^ 0.07 fm lattices. They find that the value of 
increases by about 10% with tsep, and at the same time the statistical error increases by 
a factor of 5. They use a linear extrapolation in tsep to reduce the effect of excited-state 
contamination and conclude that for the interpolating operators used tsep > 1-1 fm is 
needed to approximate the asymptotic value. In any case, one should include the excited 
states explicitly in the analysis of the matrix elements as demonstrated in Refs. [^ [S3]. 

The uncertainty in the lattice determinations of qa, which still do not fully include 



all systematic errors discussed in Sec. 5.2, is much larger than the experimental one. 



limiting its utility as a probe for physics beyond the Standard Model. Our conclusion is 
that a combination of high statistics, use of multiple tsep and investigation of correlators 
with different overlap of source with ground versus excited states will be needed to 
extract the matrix elements with high precision. A promising direction for reducing 
the statistical error in qa is to use a simultaneous chiral extrapolation of the octet 
baryons since the axial charges of the E and S baryons are calculated with significantly 
smaller errors [97]. A major limitation to testing whether excited-state contamination 
is a significant factor in the underestimate of qa is the computational resources needed 
to simulate close to (and eventually at) the physical light-quark masses, high statistics 
and extrapolations to the continuum limit. The U.S. national report on the future of 
extreme-scale computing [98] has made the high-precision calculation of qa a milestone 
to achieve, so we anticipate steady improvement in lattice estimates of all such matrix 
elements with increasing computational power. 

5.4 Nucleon induced-pseudoscalar charge Qp 

There has been renewed interest in the induced-pseudoscalar form factor gp{q^), de- 



fined in Eq. 5b, due to the recent MuCap Collaboration [99] high- precision experiment 
studying ordinary muon capture (OMC) by protons, — )■ u^n. We define the induced- 
pseudoscalar coupling as 

9p = ^Mg' = 0.88my, (37) 

where is the muon mass. Improved calculations of electroweak radiative correc- 
tions [100] allow precise extraction of the form factor from these experiments. The new 
MuCap experiment yields Qp = 7.3 ±1.1 [9911100] , which is consistent with the value pre- 
dicted by heavy-baryon chiral perturbation theory g*j^^'^ = 8.26 ± 0.16 [nJT]. However, 
it is much smaller than the earlier world average for OMC, [sflP^^Jave = 10.5 ± 1.8 given 
in Ref. [102] . and the value obtained from a TRIUMF experiment with radiative muon 
capture (RMC), fi~p t'/i^T, which gave g'^^'^ = 12 A ± 1.0 [103] . After reanalyzing 
the TRIUMF data, Clark et al. [M] found g*p = 10.6 ± 1.1. When combined with the 
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Figure 4: (Upper row) The axial charge versus from A^/ = 2[8ll|82l[83l|8ll[85] (left) 



and 



/ 



2 + 1 



[92] (right) calculations with different types of 0(a)- 
improved fermion actions. The filled symbols and solid errorbar (open symbols and 
dashed errorbar) denote results taken from published papers (the latest lattice proceed- 
ings). (Lower panel) Comparison of the published values of qa after chiral and continuum 
extrapolations with experimental measurements [96j (vertical band). The solid lines in- 
dicate statistical error while the dashed lines include systematic errors. Lin et al. [97] 
find that an SU(3)-constrained fit to the qa for octet baryons reduces the statistical error 
in the chiral extrapolation. This is illustrated by the larger errors in the LHPC result 
compared to those in Ref. [97], which are obtained using similar lattice parameters. 
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Figure 5: (Left) The induced-pseudoscalar charge of the nucleon Qp from experimental 
measurements [103^ 1102^ 1104^ [UUl IIUU] and an earher estimation from HB^PT |101j . 
(Right) Comparison of lattice estimates of Qp using the DWF fermion action [55lll05plU6] 
with MuCap data. 



new MuCap result, the world average is 8.7 ± 1.0 |100] . 

There have been few calculations of the induced charge g*p in lattice QCD. Unlike qa, 
we need to calculate the form-factor at several to extrapolate gp{q^) to the same as 
those probed in experiments. The Nj = 2 [85] and 2 + 1-flavor |106j DWF calculations 
(Tablejlj) of gp evaluated at (g^ = 0.88m^) by studying the momentum dependence of the 
axial matrix elements give 7.7(1.0) and 6.6(1.2), respectively. These central values are 
about la smaller than the world-averaged MuCap estimate, as shown on the right-hand 
side of Fig. |5] 

Direct calculations of the pseudoscalar charge gp defined in Eq. 5d have not been 
done using LQCD due to the lack of experimental motivation. One technical chal- 
lenge has been removing the contribution of the pion pole to the amputated vertex in 
the calculation of Zp in RI-MOM schemes. This has recently been overcome by using 
non-exceptional momenta in the external quark legs. We, therefore, expect to provide 
estimates for gp at the same level of precision as gT- 



5.5 Nucleon tensor charge gx 

The tensor charge gx is the zeroth moment of transversity, and can be studied through 
processes such as SIDIS (semi-inclusive deep inelastic scattering). The HERMES and 
COMPASS experiments |107l I108[ 1109] presented their first estimates of gr from data 
collected at = 2.4 GeV^. Experimentally, to estimate gr one first extracts the con- 
tribution of individual quarks as a function of the quark momentum fraction x at a 
particular Q^. To obtain the contribution of each quark, the results, estimated from mea- 
surements at a finite number of values of x, are integrated over the full range < x < 1. 
The isovector tensor charge is then given by the difference between the up and down 
quark contributions. Since this analysis requires data over the full range of x, and the 
low- a; and high-x values are not well known, improvements in precision await future ex- 
periments. Current extracted numbers are highly model-dependent. Combining SIDIS 
(HERMES and COMPASS) results with Belle e+e" analysis [M] of data collected 
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at = 110 GeV^, the best experimental estimate of gr at = 0.8 GeV^ (instead of 
= 0) is 0.771^31- 

There are also estimates from purely theoretical models. These include the Nambu- 
Jona-Lasinio model |111] and the chiral-quark soliton model |112j : unfortunately, they 
are not consistent with each other. Estimates from QCD sum rules |113j have a large 
uncertainty. 

There are several LQCD estimates of qt, and we review those listed in Table [1} The 
QCDSF collaboration's 2-flavor calculations with clover fermions [5 ^ I114j . over a large 
range of pion masses (170-1170 MeV) and 3 lattice spacings [B2], show a mild increase 
in qt with M^. RBC's 2-flavor DWF calculation |85] shows a similar trend and gave 
griMS, 2 GeV) = 0.93(6) after extrapolation to the physical pion mass and using Zt 
calculated nonperturbatively in the RI-MOM scheme. These results are summarized on 
the top of Fig. |6| 

The 2+1-flavor results from the LHPC [HE SO] and RBC/UKQCD [ST] collaborations 
are summarized in the second of Fig. [6] RBC/UKQCD used DWF for both dynamical 
and valence quarks, while LHPC used the mixed-action approach, DWF on a 2+1-flavor 
staggered (asqtad) gauge ensemble. The lattice spacings in the two calculations are 
similar, 0.1224 fm and 0.114 fm; thus, we expect similar lattice-discretization errors. 
The range of pion masses explored is also comparable, 290-760 MeV by LHPC and 330- 
670 MeV by RBC/UKQCD. Both find the dependence on the pion- mass to be small 
except at the lightest pion-mass points, 290 and 330 MeV, respectively. At these points, 
the observed downward dip could be indicative of the onset of chiral logs; however, it 
is not yet clear whether these points suffer from finite-volume and excited-state effects. 
Further studies at lighter pion masses are needed to resolve these issues. 

To extrapolate the tensor charge to the physical pion mass, we employed the heavy- 
baryon chiral perturbation theory formulation |115l I116j . The resulting formula for qt 
contains one low-energy constant and two scales at lowest order in chiral logs jll6j . Anal- 
ogous formulae for the other twist-two matrix elements, the quark momentum fraction 
(x) and helicity distribution function (Ax), which can be obtained from x-dependent 
measurements of polarized and unpolarized form factors, work well in describing the 
lattice data. We, therefore, analyzed the combined RBC/UKQCD and LHPC qt data 
using the HB^PT ansatz. The fits are highly sensitive to the data points selected, since 
the chiral log is sensitive to only one point at the lightest M,r- This lack of sensitivity 
to chiral logs is illustrated by the linear fit to the 2 + 1-fiavor data shown in Fig. |6} It 
fits all points except the one at lowest and gives qt = 1.05(2). Given this lack of 
sensitivity to the chiral-log term and the high pion masses used relative to the expected 
range of validity of this order of HB^PT, we have little reason to believe that such an 
extrapolation is well controlled. We include in Fig. |6] a HB^PT extrapolation consistent 
with the data in order to illustrate the relative size of the chiral log, which may be quite 
large and appear at pion masses not much smaller than those currently available. For 
our best estimate we use qt = 1.05(35) where the central value is from the linear fit and 
the uncertainty includes the systematic error associated with the extrapolation in M^. 
Clearly data at smaller are needed. 
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Figure 6: Summary of LQCD estimates of using Nj = 2 |114t [S21 ES] (top) and 
Nf = 2+1 [STlinn] (bottom) 0(a)-improved fermion actions. Two chiral extrapolations of 
qt are shown using the combined RBC/UKQCD DWF data on their 2.7-fm ensemble [87] 
and the LHPC mixed-action data [90]. The value at the physical pion mass from the 
linear fit is shown by the red diamond. Fits using the HB^PT ansatz are very sensitive 
to removing points at large so no error band is shown. For the data, the filled 
symbols and solid error bars (open symbols and dashed error bars) denote results taken 
from the published papers (the latest lattice proceedings). 
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5.6 Nucleon scalar charge gs 

The nucleon isovector scalar charge gs has not been analyzed in lattice calculations, in 
contrast to its isoscalar partner, the scalar density (or the nucleon-cr term). There are 
no experimental measurements of this quantity, and theoretical estimates |117j (from 
different model approximations) give rather loose bounds: 0.25 < (75 < 1. Our pre- 
liminary lattice calculations show that gs has the noisiest signal compared to the other 
matrix elements discussed previously and will, therefore, drive the size of the statistical 
ensemble required for precision studies of all the matrix elements. 

To get a first estimate of gs, we have performed calculations on two sets of gauge en- 
sembles. The first uses the anisotropic clover lattices generated by the Hadron Spectrum 
Collaboration (HSC) |118[ I119j with pion masses ranging from 390 to 780 MeV. The 
second uses Nf = 2 + 1 asqtad ensembles but calculates matrix elements with domain- 
wall valence quarks with M.^ G {350, 700} MeV. The number of configurations analyzed 
range between 200 and 650. These results are summarized in Fig. [7j The error bars 
shown are statistical. 

There is no clear guidance on how to perform a chiral extrapolation to the physical 
pion mass since the data show no evidence for chiral logs. We, therefore, made fits 
assuming a behavior linear or constant in on the full and different subsets of the 
data. In Fig. [7], we show two fits, a linear one using all the data and a constant fit to 
the five lightest values. The extrapolated value from such fits to different subsets 
of data obtained by removing the points corresponding to the heaviest and lightest 
varies between 0.6-1.0. We take the mean as the central value and 0.2 as an estimate of 
the error associated with the mass extraploation. 

In addition to the large statistical error, there is significant uncertainty in the estimate 
of the renormalization constant Zs- We have used the tadpole- improved tree-level value 
Zs = ^0, where the tadpole factor uq is the fourth root of the expectation value of the 
1x1 Wilson loop. For the HSC and DWF ensembles, Uq = 0.945 and 0.938, respectively. 
A recent nonperturbative estimate of Zs for the DWF action on lattices with a similar 
cutoff a as in our calculations, converted to the MS scheme at 2 GeV, gives 0.65 [120j . 
We expect a value closer to unity due to the smearing of links in the formulation of the 
lattice actions we use. Nevertheless, based on current nonperturbative estimates with 
different actions and link smearings, our estimate is Zs = 0.7(2) in the MS scheme at 
2 GeV. Using this value would lower gs by about 25%. We take the uncertainty in Zs 
into account by doubling the error estimate, and use gs = 0.8(4). 

5.7 Lattice estimates of tensor and scalar charges for e^^ 

LQCD calculations of gs and gx discussed in the previous sections, while theoretically 
clean, require reducing a number of systematic errors. Our current understanding is that 
finite-volume effects are small for M^rL > 4, and there is little evidence for discretization 
errors at current statistics; contributions due to excited states are smaller than statistical 
errors once the time separation tscp > 1.2 fm for the current source operators and lattice 
parameters; and chiral extrapolations gives rise to the biggest uncertainty in the current 
data as shown in Figs. |6] and [7j Thus, we need high-statistics calculations on large 
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Figure 7: The scalar charge gs from Nf = 2 + 1 anisotropic clover and DWF/asqtad 
lattices. We also show a linear fit to the full data set and a constant fit to the data at 
the five smallest values of M^. The extrapolated values are shown using red diamonds. 
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lattices with light-quark masses close to the physical values. Lastly, nonperturbative 
calculations of renormalization constants are essential. 

Based on the above analysis, preliminary LQCD estimates are 

^t(MS, /i=2 GeV) = 1.05(35) , ^^(MS, /i=2 GeV) = 0.8(4). (38) 

These are used in the next section to explore bounds on new physics at the TeV scale. 
We emphasize that our focus at this point, given the preliminary nature of the estimates, 
is on the variation in the bounds under different scenarios of reduction of errors in lattice 
calculations. 



6 Impact of lattice results on phenomenology 



In Section |4.2| while studying the low-energy phenomenology of ^s^t, we ignored the 
uncertainty in the charges gs,T- Clearly, the impact on es^r of future 10~^-level neutron 
measurements of b, b^, and b^ — b depends on how well we know the nucleon matrix 
elements gs,T- Since gs,T always multiply factors of the short-distance couplings in 
physical amplitudes, they determine the slope of the bands on the e^-er plane represented 
in Figs. |2]and[3j Moreover, if one accounts for the uncertainty in qs^t the bands Figs. |2] 
and |3] acquire additional theory- induced thickness and their boundaries are mapped into 
characteristic "bow-tie" shapes. We illustrate this in Fig. [8} assuming experimental 
sensitivities in b and b^ — b at the 10~^ level. For the scalar and tensor charges we use 
in the left panel the ranges quoted in Ref. [20] (based on earlier quark- model estimates): 
0.25 < Qs < 1-0, 0.6 < qt < 2.3; while in the right panel we use the lattice estimates 
Qs = 0.8(4) and qt = 1.05(35), corresponding to Sgs/gs ~ 50% and Sgx/gr ~ 35%. 
Comparing these plots to the ones in Fig. [2] the loss of constraining power is quite 
evident. Especially in the left panel one sees that the impact of neutron measurements 
is greatly diluted. 

In Fig. |9] we summarize the low-energy constraints on es,T, taking into account the 
effects of hadronic uncertainties. We plot the combined 90% C.L. regions in the e^- 

plane allowed by the current limit on 69+ and future 10~'^-level measurements of b 
and b^ — b in neutron decay. The different curves reflect four different scenarios for 
the hadronic matrix elements: the outer-most curve corresponds the range of Ref. [20j . 
while the three inner curves correspond to lattice results with current central values from 



Eq. (38) and three different uncertainties: Sgs/gs ^ {50%, 20%, 10%} with dgx/gx = 
2/3 Sgs/gs (this choice assumes that the ratio of fractional uncertainties in gs and gx 
will remain approximately constant as these uncertainties decrease). 

The confidence intervals on es,T are obtained using the so-called R-Fit method, as 
described in Ref. |121] . In this approach the QCD parameters gs,T are bound to remain 
within allowed ranges determined by the lattice calculations and estimates of systematic 
uncertainties (in the case at hand the ranges are 0.4 < gs "£ 1-2 and 0.7 < gr ^ lA). 
The chi-squared function 

X (e5, er, gs, 9t) = 2_^[ ^ 1 (39) 
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is then minimized with respect to gs,T (varying qs^t in their allowed ranges), leading to 

X^(e5, ct) = min^g J, ct, gs, 9t) ■ (40) 

Finally, the confidence intervals on es,T are deduced applying the standard procedure [96] 
to x^{es, ex), with an effective number of degrees of freedom given by min(iVobs — -^V^, N^), 
where A'obs is the number of experimental constraints, A'g = 2 is the number of QCD 
parameters {gs,T), and A^^^ = 2 is the number of parameters we wish to constrain {es^r)- 
From Fig. |9] several clear messages emerge: 

• Hadronic uncertainties in qs^t strongly dilute the significance of new 10~^-level 
experiments. Experimental progress without theoretical progress will not lead to 
competitive constraints on the short-distance scalar and tensor interactions. 

• Our preliminary lattice results (curve labeled by Sgs/gs = 50%) already provide a 
significant improvement over previous knowledge of gs,T summarized in Ref. |20j . 

• In order to fully exploit the constraining power of planned 10~^ measurements of 
b and b^, the uncertainty on gs should be reduced to 20%. Improvement beyond 
this level would not significantly increase the constraining power (see difference 
between the curves labeled as Sgs/gs = 20% and Sgs/gs = 10%). 
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Figure 8: Left panel: 90% C.L. allowed regions in the c^-ct plane implied by (i) the 
existing bound on 60+ (green horizontal band); (ii) projected measurements of b and b^ — b 
in neutron decay (red and blue bow-tie shapes) at the 10~^ level; (iii) hadronic matrix 
elements taken in the ranges 0.25 < gs < 1-0, 0.6 < gx < 2.3 [20]. Right panel: same as 
left panel but with scalar and tensor charges taken from lattice QCD: gs = 0.8(4) and 
gr = 1.05(35). The effective couplings es,T are defined in the MS scheme at 2 GeV. 
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Figure 9: Combined 90% C.L. allowed regions in the es-er plane based on: (i) existing 
limit on 6o+ from O"*" — )■ O"*" nuclear decays; (ii) future neutron decay measurements 
with projected sensitivity of 10~^ in b and b^, — b. The four curves correspond to four 
different scenarios for the hadronic matrix elements: 0.25 < gs < 1-0, 0.6 < Qt < 2.3 
as quoted in Ref. [20]; lattice results with current central values from Eq. (38) and 
Sgs/gs = 50%, 20%, 10% with Sgr/gr = 2/3 6gs/gs (this choice assumes that the ratio 
of fractional uncertainties in gs and gr will remain approximately constant as these 
uncertainties decrease). The effective couplings es,T are defined in the MS scheme at 
2 GeV. 



7 Collider limits 

The contact interactions probed at low energy can also be directly probed at high-energy 
colliders. The rate, however, depends on whether the particles that generate the 4-fermi 



interaction are kinematically accessible at the collider energies. We begin in Section 7.1 



under the assumption that the scalar and tensor interactions remain point-like at TeV 



scale energies. Then in Section 7.2 we derive a relation between eg and the production 



cross-section, Eq. 54 , when the scalar interaction is generated by the exchange of a 



resonance that is kinematically accessible at the LHC 



7.1 Model-independent limits 

Assuming that the scalar and tensor interactions remain point-like at TeV-scale en- 
ergies, we can employ the operator formalism to put bounds on es,T,p from collider 
physics. SU (2) gauge invariance implies that es,T,p control not only charged- currrent 
processes but also the corresponding neutral-current versions, as the weak-scale effective 
Lagrangian includes terms proportional to (eg — ep)eiieLdLdfi, (e^ -|- ep)eReLUjiUL, and 
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^T(iRcr^^eLUjia^uUL- Exploiting this property, from an early CDF analysis [122J of contact 
interactions in pp — )■ e'^e' + X, after matching the different conventions for the effec- 
tive couplings, we obtain the 90% CL. limit \es\ < 0.135. There are a number of LHC 
searches for contact interactions, specifically in dijet |123[ 11241 1125] and dimuon \l'2b\ 
final states. All of these studies, however, focus only on specific vector-like interactions 
and do not consider scalar (i.e, helicity flipping) contact interactions. 

Here we focus for definiteness on the charged-current part of the scalar and tensor 
effective operators. These contact interactions fall into the signature class of collider 
searches for an exotic W gauge boson, since they both can contribute to the signature 

— )■ ez/ + X. We will use the analyses and results of searches for this process to obtain 
bounds on and ey- In the limit m; = the analysis is simplified, since these operators 
do not interfere with SM processes. We do include the interference between the scalar 
and tensor interactions, which does not vanish in the chiral limit. The relevant part of 
the effective Lagrangian is given by 

C = -^KdM(ePz.^^e)-^Kd(na^'^PLrf)(ea^.PL^^e) (41) 

where a^'^ = ^[7'^, 7'^]/2, and f]s,VT = =t denotes the sign of the coefficients of the 
scalar and tensor operators. The relations between and the effective couplings 
es,T at yU = 1 TeV are given by = 2r7st>^/A| and er = rixv"^ /K\. Note that since 
collider searches set limits on the effective couplings es^T at the high renormalization scale 
/i = 1 TeV, a direct comparison with the low-energy constraints requires an appropriate 
rescaling down to the hadronic scale. Using the one-loop anomalous dimensions for scalar 
and tensor operators (see |127j and references therein), the one-loop beta function for 
the strong coupling constant, and including the appropriate heavy quark thresholds, we 
find in the MS scheme es{l TeV)/es(2 GeV) = 0.56 and €7(1 TeV)/eT(2 GeV) = 1.21. 
We will use these factors to rescale the collider limits and compare them to low-energy 



limits in Figs. 10 and 11 



To determine the transverse mass distribution of the electron-neutrino pair we start 
with [HH] 

= -^yl^h^nu?) (42) 

dy dy' dm% QA-ns^ ^ Xi X2 

tj 

where i and j are summed over the initial partons (with parton distribution func- 
tions (PDF) fij and momentum fractions Xi^2)i and y,y' are the rapidities of the elec- 
tron and neutrino. One finds xi = mT{ey + e'^)/2y/s, X2 = rriT^e'^ + e~'^)/2y/s. 
We also used the observation that the transverse mass of the electron and neutrino, 



rriT = \/2E^E!^(l — cos ^(peu) (where E^^ is the transverse energy of the electron or 
neutrino, and ^(pev is the azimuthal angle between the two leptons), is simply uit = 2pT 
at leading order, where px is the transverse momentum of the electron. 

To leading order (LO), the contributions of the color- and spin-averaged scalar and 
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tensor matrix elements to pp eu + X, including the interference term, is 



{p-p'){k-k') 



l^^[ip-kW-k')-{p.k'W-k)] 



2 1 
3A| 
16 1 

+ y ^ [2(P ■ k){p' ■ k') + 2(p • k'){p' ■k)-{p- p'){k ■ k')] 



(43) 



where p,p' are the momenta of the incoming partons, and k, k' are the momenta of the 
electron and neutrino. The interference term is antisymmetric under k < — > k', so it does 
not contribute to the transverse mass distribution obtained by integrating over y and y'. 
After some substitutions this expression becomes 



{\M\') 



6A| 



VsVt m^s , 4 ^2 



1 - 



rrin 



(44) 



with s = X1X2S and rriT is the transverse mass of the lepton-neutrino pair. 
Next we need the cross-section with my greater than a threshold m^^cut- 



Using (42) 



and (44), one finds 



airriT > mx. 
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(45) 



where the sum over i,j = u,d and i,j = u,d has been done, r = 0:1X2, and yp = 
0.5 ln(a;i/x2). 

The ATLAS and CMS experiments have searched for new physics in pp ^ eu + X 
by looking for an excess of events at large transverse mass |129l I13U] . The CMS study 
analyzes 1.03 fb~^ of data for the electron final state, and we begin by following their 
analysis in setting limits on the scalar and tensor interactions. The CMS search window 
is defined by specifying a cut on rriT and counting the number of events detected with 
transverse mass larger than the cut. Specifically, they looked for the production of a 
heavy W with decay W — )■ eu by searching for events having transverse mass above a 
variable threshold mT,cut; finding 1 event for m^cut = 1 TeV and 1 event for m^-cut = 
1.1 TeV. In general the limit on the number of expected signal events depends on the 
expected background, nb, which for this search is quoted to be nb = 2.2 ± 1.1 events for 
"^T,cut = 1 TeV and nb = 1.4 ± 0.80 events for mT,cut = 1-1 TeVp^ 

To set a limit we follow Ref. |130j and use Bayesian statistics with a flat prior 
in the signal ng. The likelihood function L{n\ns) is given by the Poisson distribution 
for n detected events with ng signal and rib background events expected. The expected 
number of signal events is given by Ug = eaC, where a is given by (45 ), £ is the integrated 
luminosity, and e is the detection efficiency times the geometric acceptance. Ref. |130] 



These rib values are taken from Table 1 of Ref. 

1.15 for mT,cut 



Different central values for rib appear in 
1 TeV. The two sets of nb are consistent 



Figure 2 of Ref. |130j : for example rih 
within the quoted error bars and lead to minor (5%) differences in the bounds on es,T 
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quotes the signal efficiency for a W to be 80%. Their earher analysis (Ref. |131j ). 
based on 36 pb~^ of data, quotes the product of the geometric acceptance and detection 
efficiency as being greater than 64% in the W mass range of interest. In the absence of 
a detector simulation for our signal, in what follows we will assume our signal has a 80% 
detection times geometric acceptance efficiency. 

The credibility level 1 — a for a flat prior in the signal is then derived from [96] 



fn"^ dus L(n\ns) , ^ 

which is equivalent to [96] 

a = e-^- ^^n.=oM^up + n^r ^^^^ 

Z^m=0 m! b 

To set a limit on Sup, we choose the lower value of mT,cut = 1 TeV in order to maximize 
the signal rate. Then for nb = 2.2 expected background events and n = 1 event detected, 
one finds that s^p = 3.0 at the 90% credibility level. Dividing by e, we obtain a 90% 
upper credibility limit of 3.7 produced signal events. 



In Fig. 10 we show the corresponding limits on and (red, solid curve), using Eq. 
(|45[ with 1.03 fb-i of integrated luminosity at = 7 TeV. LO MSTW 2008 [I32] PDFs 
are used and evaluated at = 1 TeV^. We also checked that limits obtained using the 
CTEQ6 PDF set |133j are quantitatively in good agreement. When only one of these 
operators is present, the bounds on es^x correspond to As > 2.5 TeV and > 2.7 TeV. 



As illustrated by Fig. [10} our LHC bound on is a factor of 7 stronger than the old 
CDF limit, but about 4 times weaker than the bound from nuclear beta decay. Similarly, 
our new collider bound on is a factor of 3 weaker than the bound from radiative pion 
decay. 

We have performed a parallel analysis using the ATLAS results [129j on W search. 
Use of the ATLAS results requires an extra step, since their quoted efficiency for a 
given rriT^cut includes the fraction of total W — )■ eu events with rriT > "^T.cut- After 
determining this fraction with a leading-order calculation, we infer the ATLAS detection 
times geometric acceptance efficiency for m^cut = 1 TeV to be 80%, the same as quoted 
by CMS for their experiment. Using then the fact that for ttlt > 1 TeV ATLAS observes 
n = 1 event with an expected number of background events rib = 0.89(20), we find that 
the ATLAS limits on es,T differ from the CMS ones only at the 5% level, well within the 
uncertainties of our leading order calculation. We also estimate that the bounds on es,T 
can be reduced by at least a factor of 2 once the full data set collected at y/s = 7 TeV 



is analyzed (see dotted, gold line in Fig. 10). We expect that stronger limits can be 



obtained by a combined analysis of the ATLAS and CMS data, which goes beyond the 
scope of this work. 

To obtain projected limits at higher luminosities and ^/s = 7 TeV or ^/s = 14 TeV, 
we repeat the same LO analysis, assuming the same 80% detection times acceptance 
efficiency for the signal as before. We choose an aggressive cut to make the expected 
background small. The location of the cut on the transverse mass m^cut ~ TeV is 
estimated by computing at tree-level the transverse mass distribution of the dominant 
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SM physics background, due to the production of a high-p-r lepton from an off-shell W, 
and finding the value above which the expected background is less than 1 event. At 
y/s = 7 TeV and with an integrated luminosity of 10 fb~^, we find that the number 
of background events drops below one for m^cut = 1-5 TeV. At y/s = 14 TeV, with 
'^T.cut = 2.5 TeV and an integrated luminosity of 10 fb~^ we find 0.5 background events 
expected above the cut. We also consider an ultimate luminosity of 300 fb~^, finding 
that for mj' > 4 TeV there are 0.3 expected background events. We therefore impose 
"^T,cut = 2.5 TeV (4 TeV), and assume an integrated luminosity of 10 fb~^ (300 fb~^). 
To set a limit we will assume that no events are found, which is consistent with less than 
1 background event expected. From Eqs. (46) and (47) we then obtain a 90% Credibility 
Limit of 3 produced events. The anticipated joint 90% credibility level limits on and 
St from LHC at ^/s = 14 TeV are shown in Fig. 
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From an inspection of Fig. 11 we find that at high luminosity and center-of-mass 



energy the expected improvement in the limits are nearly an order of magnitude com- 
pared to the existing collider limits. Even with only 10 fb~^ taken at 14 TeV we expect 
the limits to improve substantially from the current collider limit. At these energies and 
luminosities the bound on es from the LHC will become stronger than anticipated future 
bounds from low-energy experiments. This conclusion is illustrated in Fig. [TT| where we 
also overlay the projected low-energy bounds presented in Fig. |9| 

We expect these projected limits can be tightened, since we have chosen hard cuts to 
reduce the expected leading background to below one event, at the cost of significantly 
cutting into the signal. Optimizing the choice of the cut to maximize the sensitivity to 
the contact interactions will require including additional backgrounds, such as QCD and 
top quarks, and more generally, a better understanding of the systematic errors involved. 

Our analyses can certainly be improved. Our estimate of the detection times accep- 
tance efficiency was borrowed from the estimate for a W signal from |130] . Obviously 
a detector simulation of the signal will provide a better estimate of this factor. The 
theoretical error on the signal rate will be further reduced once the next-leading-order 
QCD contributions are known. 

Finally, the interaction Lagrangian (41) can be generalized to include interactions of 
the electron with neutrinos of all flavors, with a corresponding generalization of es,T 
eg J' where a G {e,/i, r}. Because the final states with different neutrino flavors do 
not interfere and neither do the scalar and tensor interactions after integrating over the 
rapidity distributions, the derived and projected bounds shown in Figs. 10 and 11 now 

apply to the quantitiesA/^^(egy)^. 



7.2 Scalar resonance 

A larger signal rate is obtained if the particle that generates the scalar interaction is 
kinematically accessible at the LHC. In this case there can be a direct relationship 
between and the production cross-section and mass of the resonance, as we now 
demonstrate. 

We assume that after electroweak symmetry breaking there is a charged scalar 0+ of 
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Figure 10: Joint 90% CL limit on and implied by: (i) cmrent bounds from nuclear j3 
decay 0+ — )■ 0"*" and radiative pion decay (blue, dashed); (ii) CMS search [130j in the channel 
pp ^ ev + X ai y/s = 7 TeV with 1.03 fb^^ of data. The limit is obtained by requiring less 
than 3.7 ez/-produced events having rriT > 1 TeV (red, solid). LO MSTW 2008 |132j parton 
distribution functions are used; (iii) projected LHC searches in the channel pp ^ ev + X at 
y/s = 7 TeV with 10 fb~^ of data (gold, dotted). The limit is obtained by requiring less than 3 
ei/-produced signal events with > 1.5 TeV and assuming that no events are observed. The 
cut is chosen to reduce the expected leading background to be below 1 event. The effective 
couplings es^T are defined in the MS scheme at 2 GeV. 
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Figure 11: Projected joint 90% credibility level limit on es and from the LHC at ^/s = 14 
TeV, obtained from requiring less than 3 ez^-produced signal events with: (i) uit > 2.5 TeV and 
10 fb~^ of integrated luminosity (solid, red ellipse); and (ii) uit > 4 TeV and 300 fb^^ (dashed, 
yellow ellipse). Cuts are chosen to reduce the expected leading background to be below 1 event. 
To obtain the projection it is assumed no events are found. Same PDFs are used as in Fig. 



10 Note the change in scale between these two figures. Anticipated bounds from low-energy 
experiments and reduced LQCD uncertainties, redrawn from Fig. [9j are shown for comparison. 
The effective couplings es t are defined in the MS scheme at 2 GeV. 
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mass m, with the following couplings to first-generation quarks and leptons: 



£ = XsVud(f)'^ud + XpVud(j)'^u^5d + Xi(j) ePiVe + h.c. , (48) 
where 0~ = (0^)*. At low energies the exchange of 0"'" generates a scalar operator with 



es = 2A5AZ— (49) 



and a pseudoscalar operator with 



m 



ep = 2XpXi— . (50) 



To proceed, at leading order the cross-section for the on-shell production of (p, which 
then decays to lu (of a given sign), is given in the narrow-width approximation by 

a-BR = Af(A| + A?.)|K/^^(r) (51) 

with r = m^/s, L(r) = dxfg{x)fg{T/x)/x, and where is the total decay width of 
(p. Next, note that since may decay to other particles (not just to lu and ud) , 

r<^>r, + r, = (A2 + 2iv,(A| + A|)|KdP) ^ (52) 

with Nc = 3. Next note that m/T^ < 167r/(Af + 2A/'c(A| + A^lKdP), and then use 
the arithmetic-geometric inequality A/2iVcA(^A; < |(Af + 2Nc{X'g + Xp)\Vud\'^), where 
A^ = ^/X^~+~X^\Vud\, to finally obtain our main result of this subsection, 

^■■^■^ ^ ^7fc(^)^^<^> '''' 

Because of the severe constraint imposed by vr — )■ eu, the coupling Ap of to the 
pseudoscalar quark scalar density must be significantly suppressed. In the limit ep <^ 
one then has 

a-BR < ^-^|es|rL(r) (54) 

This expression can be rearranged to obtain a lower bound on es, that is stronger after 
summing in L over both charged-particle final states. The bound depends only on r and 



a ■ BR. Figure 12 shows the bound as a function of r for several choices of a ■ BR that 



will be probed by the LHC Equivalent limits are shown in Figs. 13 and 14, where we 
show the dependence of the bound on the mass of the resonance, for several values of 
a ■ BR and at ^/s = 7 and 14 TeV. LO CTEQ6 |133] parton distribution functions are 
used for all these figures. 

We have shown that if a signal is observed in pp — )■ e+missing energy (MET)+X, then 
a lower bound on es can be obtained, provided the signal is due to the on-shell production 
of a scalar, which decays to an electron and missing energy provided by an electron 
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Figure 12: Projected lower bound on |e5| (at /i = 2 GeV), for a discovery cross-section 
of a{pp ^ e+ MET + X) = 10 fb (blue, solid), 1 fb (red, dashed) and 0.1 fb (black, 
dotted), as a function of r = m? / s. Shaded region (green) shows the current experimental 
exclusion on from 0+ — )■ 0+ nuclear (3 decay. 




m(TeV) 

Figure 13: Projected lower bound on |es| (at = 2 GeV) for y/s = 7 TeV and a discovery 
cross-section of a{pp e + MET + X) = 5 fb (blue, solid), 1 fb (red, dashed) and 0.5 
fb (black, dotted). Shaded region (green) shows the current experimental exclusion on 
from 0^ — t- O"*" nuclear /3 decay. The bound scales linearly with cr ■ BR. 
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neutrino, and whose pseudoscalar couphng to quarks is suppressed compared to its scalar 
coupling. That the resonance couples to an electron neutrino is important in deriving 
the above relation to es, since at linear order in the e's, neutron-decay experiments do 
not probe couplings to other neutrino flavors e^^^^'^^. 

Further confidence that the signal is due to the production of an on-shell particle 
can be established by the detection of an edge in the transverse electron-neutrino mass 
distribution, and through the detection of a resonance in dijets. The only additional 



theoretical assumption used to obtain the lower limit (54) is that the charged resonance 



is interpreted as a scalar and not as a vector or a tensor. Measurements of the rapidity 
distribution of the electron should determine the spin of (p. 

If a signal is discovered in pp — > e+MET+X then neutron and nuclear /3-experiments 
will be crucial in order to pin down the properties of the resonance and of the MET. 
As an illustration, suppose the measured cross-section and mass imply an €5 in excess 
of existing neutron decay bounds, then either: (i) the resonance does not have spin 0; 
or (ii) as already described above, since the outgoing neutrino flavor is not identified, 
the relationship between the cross-section and can be undone simply if the scalar 
couples preferentially to muon and tau neutrinos rather than to the electron neutrino; 
or (iii) there are additional scalars at the TeV scale or a scalar contact interaction, such 
that partial cancellations occur in summing the multiple contributions to e^. 



8 Discussion 

It is anticipated that the next generation of neutron /3-decay experiments will increase 
their sensitivity to BSM scalar and tensor interactions by an order of magnitude, through 
improved measurements of the neutrino asymmetry parameter B and the Fierz inter- 
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ference term b (see Figs. |2| |3| and [9]). In order to assess the impact of these future 
experiments, we have performed a comprehensive analysis of constraints on scalar and 
tensor BSM interactions from a broad range of low-energy probes (neutron decay, nuclear 
decays, pion decays) as well as collider searches. 

Extracting bounds on scalar and tensor BSM couplings from neutron and nuclear 
beta decays requires knowledge of the nucleon scalar and tensor form factors at zero 
momentum transfer. In this paper we have provided the first lattice-QCD estimate 
of the scalar form-factor, gs = 0.8(4), and a new average of existing tensor form-factor 
results, gx = 1.05(35). We find that to fully exploit the increased experimental sensitivity 
will require understanding the lattice-QCD estimates of the proton-to-neutron matrix 
elements at the level of 10-20% (see Fig. |9|. To do that will require analyzing a few 
thousand samples at each value of the simulation parameters using a combination of 
decorrelated lattices and multiple source points on each lattice, improvements in source 
and sink interpolating operators for nucleons, and simulations close to physical light- 
quark masses. With the anticipated increase in computing power and resources, we 
estimate calculations will reach this precision in 2-4 years. 

In our survey of probes of BSM scalar and tensor interactions, we have found that 
the currently strongest bounds arise from nuclear (3 decay (65) and radiative pion decay 
(ey), probing effective scales > 4.7 TeV and > 5 TeV, respectively. We also 
find that within a specific model for the lepton flavor structure of the scalar and tensor 
interactions, significantly stronger bounds arise from tt — )■ ei/ decay, a conclusion in 
agreement with previous literature. 

We have used LHC data to obtain constraints on the scalar and tensor interactions, 
finding bounds within sight of current limits obtained from low-energy measurements 
(see Fig. 10). We have also provided a preliminary estimate of expected future bounds 
from the LHC, finding that an order of magnitude improvement should ultimately be 
achievable and that the future collider constraints (associated with effective scales As^t ^ 
7 TeV) will compete with improved neutron-decay constraints based on experimental 
sensitivities 5h,5hy ~ 10^^ (see Fig. 11). Finally, if a charged resonance decaying to 
an electron plus missing energy is discovered at the LHC, we have shown how, with 
some theoretical assumptions, the production cross-section provides a lower bound on 
the scalar interaction probed at low energy (see Figs. 12, 13, and 14). 

Our analysis shows that in order to compete with upcoming collider bounds on scalar 
and tensor interactions, future neutron-decay experiments should aim at the very least 
to sensitivities 56, Sh^ ~ 10~^ in the Fierz interference term and neutrino asymmetry. 
Moreover, experiments aiming for 56, 5hv ~ 10~^ would provide an unmatched discovery 
potential for new scalar and tensor interactions, and therefore should be vigorously 
pursued. 
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A SU{2) X [/(I) invariant operators contributing to 
charged-current processes 

The building blocks to construct gauge-invariant local operators are the gauge fields 
G^, Wj^, B^, corresponding to SU{3) x SU{2)l x U{1)y, the five fermionic gauge mul- 
tiplets, 



the Higgs doublet ip 



and the covariant derivative 



^-ilo ] , (56) 



D.^Id,- igs^Gt - ig^K - . (57) 

In the above expression A"^ are the SU (3) Cell-Mann matrices, o"" are the SU (2) Pauli 
matrices, gs,g,g' sue the gauge couplings and Y is the hypercharge of a given multiplet. 

The minimal set of operators contributing to low-energy charged current processes 
can be divided into two groups: four-fermion operators 

Ojf = (/7^a"/)(g7^a'^g) (58a) 

Ogde - ile)(dq)+h.c. (58b) 
Oig = {lae)e''\qbu) +h.c. (58c) 
Ol = {laa^''e)e''\qba^,u) + h.c. (58d) 



and vertex corrections 



O^^ = i((^^eD^v^)(M7^rf) +h.c. , (59a) 
OjfJ = t{ip^D'^ay){q'y^a^q)+h.c.. (59b) 
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Moreover, the extraction of the Fermi constant from muon decay (needed for weak uni- 
versahty tests) is affected by two more operators (four-lepton and lepton-gauge vertex 
correction) : 



0(3) 

o 



II 

(3) 



2(V9tD^(TV)(l7^(T"/)+h. 



(60a) 
(60b) 



In terms of the coefficients of the above operators, the low-energy effective couphngs 
appearing in Cqc (see Eq. [2]) are given by 



V.J- 




Uij 


= 2 V, [a^f 


-\-2V- 

U 


Vr 




Uij 






Vr 




Uij 


= — ['^iq]i£ji 




Vr 


[sr 


Uij 


Vim \oiqde 


* 

Ujm 


V^ 




Uij 







.(3)1 



a 



(3) 
Iq 



Umj 



(61a) 

(61b) 
(61c) 
(61d) 
(61e) 



B Details of neutron decay distribution 

The effective Fierz interference term h and effective energy-dependent correlation coeffi- 
cients a{Ee)^ A{Ee), B{Ee) and C(^aa,aA,aB){Ee) introduccd in Eq. plare [36l |37] : 



b = foSM^^BSM (g2a) 

-aiE^) = («Lo(A) + c('^)+^^^) (l + ^en^^^e)) (62b) 

A{Ee) = (^Lo(A) + 4^) + cS^)^)(l + ^5i^)(i^e)) (62c) 

B{E,) = 5lo(A) + 0^ + 4^)^ + 1^ (6^^ + (62d) 

C^aa){E.) = c^r^^ (62e) 

CiaA){E^) = c^^^ (62f) 

Ciasm = (cr^+cr^) ■ (62g) 

In these expressions the subscript LO indicates the well-known leading-order contribu- 
tions that survive if we neglect the radiative corrections, recoil effects and new-physics 
contribution^^ 

(\\ 1-^' 4 (^^ 2A(1-A) 2A(1 + A) 

aLo(A) = ^^, ^Lo(A) = ^^^, 5lo(A) = ^^^. (63) 
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In that limit, of course A — > A in olo; ^lo and -Blo- 
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As discussed in the main text, the hnear new-physics effect due to the eR couphng has 
been included in the replacement A — )■ A = A(l — 2e/j). The only other linear BSM effects 
in the differential distribution are b^^^ and b^^^, whose expressions are shown in the 



main text, Eqs. (11). 

Radiative corrections are encoded in the function 6a {E^) |3B], while recoil corrections 
are encoded in the coefficients (.^^^^^^"■"■^"■'^'"■^ ^ ^sm ^^^^ ^sm^ whose explicit expressions 
are [37] 

_ 2A(A + /iv) Eo 
'° - " 1 + 3A2 ^^^^ 
3 + 4A/iv + 9A2 
= 1 + 3A^ 

(a) _ 2A(A + Eq 

~ 1 + 3A2 ^^^^ 



4A(3A + /xy) 
1 + 3A2 



(A) ^ /xy(l-3A) + A(7-5A) 
1 + 3A2 

^(B) _ 2A(A + /iy) 



1 + 3A2 Mn 
/iv(l + 3A) + A(5 + 7A) 



1 + 3A 



2 



(67) 



(A) _ (A - l)(A + /iv) ^0 , ^ 

" 1 + 3A2 ^^^^ 



(69) 
(70) 
(71) 



_ (l + A)(A + /iv) ^0 

" 1 + 3A2 M^. ^^^^ 

{aB) 

(/iy + 7A)(l + A) 

" 1 + 3A2 ^^^^ 

>a) ^ 3(1 - A^) 

^ 1 + 3A2 

{aA) _ (A - l)(/iy + 5A) 

~ 1 + 3A2 



(74) 
(75) 



,SM ^ me l + 2/xvA + A2 

1 + 3A2 ^ 

,SM ^ me (l + A)(/iy + A) 

1 + 3A2 ^ ' 

In the above relations represents the difference between the proton and neutron 
magnetic moments. Numerically, one has = —1.35(1) x 10^^ and b^ = —1.27(1) x 

lo-l 
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